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FOREWORD 


This  report  covers  research  performed  under  U.  S.  Army  Research  Office 
Contract  DAAG29-83-C-0001  during  the  period  October  10,  1982  to  December  20, 
1985.  The  high  pressure,  high  temperature  CARS  spectra  of  nitric  oxide  and 
water  vapor  were  Investigated  for  the  purpose  of  permitting  accurate  temper¬ 
ature  and  concentration  measurements  of  these  species  In  combustion  media. 

This  was  accomplished  by  the  acquisition  of  high  quality  CARS  spectra  under 
controlled  laboratory  conditions  of  temperature  and  pressure  which  were  used 
to  refine  and  verify  a  theoretical  model  that  predicts  the  spectra  through 
computer  simulation.  CARS  measurements  In  high  pressure  flames  were  Inves¬ 
tigated  through  the  acquisition  and  analysis  of  CARS  spectra  from  a  premixed 
flame,  operated  at  pressures  up  to  15  atmospheres. 

The  personnel  Involved  with  this  contract  have  had  responsibilities  for 
different  phases  of  the  work.  John  A.  Shirley  performed  the  experimental 
Investigations  Into  H2O  CARS  spectra  and  was  responsible  for  their  compar¬ 
ison  with  the  theoretical  models.  Similarly,  John  H.  Stuff lebeam  did  the  CARS 
studies  of  NO  and  was  also  responsible  for  the  high  pressure  flame  measure¬ 
ments,  Including  the  analysis  of  different  laser  convolution  predictions  for 
the  CARS  spectra.  Gregory  M.  Dobbs  provided  computer  programs  based  on  the 
Kataoka-Teets  laser  convolution  used  for  comparison  with  the  experimental 
data.  Alan  C.  Eckbreth  provided  the  overall  direction  of  the  experimental  and 
theoretical  aspects  of  this  work  and.  In  addition,  supplied  specific  solutions 
to  experimental  problems  that  were  encountered.  This  contract  was  monitored 
by  Dr.  Robert  W.  Shaw  of  the  Chemical  and  Biological  Sciences  Division,  U.  S. 
Army  Research  Office. 
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SUMMARY 


Under  Contract  DAAG29-83-C-0001  sponsored  by  the  U.S.  Army  Research 
Office,  the  United  Technologies  Research  Center  has  conducted  spectroscopic 
investigations  into  the  coherent  anti-Stokes  Raman  spectroscopy  (CARS)  spectra 
of  NO  and  H2O  and  investigated  the  capability  of  CARS  for  diagnostics  of 
high  pressure  flames.  CARS  is  a  remote  laser  diagnostic  technique  for 
temperature  and  species  measurements  in  hostile  combustion  environments.  As 
such  it  possesses  considerable  relevance  to  the  Army  in  the  general  areas  of 
ballistics  and  propulsion.  This  final  report  describes  the  results  of 
investigations  conducted  into  the  effects  of  high  pressure,  specifically  the 
phenomenon  of  colllslonal  narrowing,  on  CARS  spectra  from  which  temperature 
and  density  information  derive.  Experimental  studies  of  CARS  spectra  were 
conducted  in  NO  in  a  heated,  high  pressure  cell  over  the  range  300  K-800  K, 
1-33  atmospheres,  in  H2O  from  500  K-1000  K,  5-30  atmospheres  and  in  a  high 
pressure,  premixed  CO/alr  flame  from  1-15  atmospheres.  The  CARS  spectra 
obtained  were  used  to  evaluate  the  UTRC  CARS  computer  code.  In  addition,  a 
new  technique,  termed  dual  broadband  CARS,  was  developed  which  permits  the 
simultaneous  detection  of  a  multiplicity  of  species.  This  approach  removes 
the  nominal  limitation  on  CARS  of  interrogating  just  a  single  constituent  at  a 
time.  This  technique  should  prove  Invaluable  in  the  transient  environments  of 
combustion  processes  relevant  to  the  Army. 
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INTRODUCTION 


High  pressure  cumbustlon  Is  extremely  important  In  a  variety  of  practical 
applications  of  relevance  to  the  Army  such  as  propulsion  and  ballistics. 
Experimental  diagnostics  of  high  pressure  combustion  phenomena  are  Important 
to  gain  the  understanding  required  to  Improve  and  control  these  processes. 
Optical  techniques  appear  ideally  suited  to  diagnosing  such  phenomena  because 
they  are  nonintrusive ,  can  be  spatially  and  temporally  precise  and  are  not 
confronted  with  survival  considerations.  Many  however  are  inappropriate  due 
to  intolerance  to  interferences.  One  optical  technique,  coherent  anti-Stokes 
Raman  spectroscopy  or  CARS,  shows  considerable  promise  for  diagnosing  high 
pressure  combustion  phenomena.  CARS  is  spatially  and  temporally  precise  and 
its  coherent  or  beamlike  nature  is  quite  amenable  to  the  limited  optical 
apertures  typical  of  high  pressure  combustion  chambers.  Furthermore,  CARS 
exhibits  a  nonlinear  dependence  on  molecular  number  density  and  increases 
rapidly  in  signal  strength  as  the  gas  density  is  elevated.  The  intensity  of 
CARS,  however,  is  quite  llnewidth  dependent.  Collision  processes  alter  the 
spectra  by  perturbing  the  rovibronlc  llnewidth.  As  the  gas  density  Increases, 
a  variety  of  llnewidth  phenomena  occur  such  as  pressure  broadening,  shifting, 
motional  or  Dlcke  narrowing,  and  collisional  narrowing  which  complicate  the 
calculation  of  the  CARS  spectrum.  For  CARS,  the  collisional  narrowing 
phenomenon  is  Important  for  pressures  above  one  atmosphere.  Collisional 
narrowing  occurs  when  adjacent  Raman  transitions  have  been  pressure  broadened 
to  the  extent  that  they  overlap.  The  overlap  allows  communication  between 
lines  and  they  subsequently  coalesce  or  collapse  to  a  narrower  bandwidth. 
Accurate  modeling  of  high  pressure,  high  temperature  CARS  spectra  is  necessary 
to  extract  temperature  and  density  information  and  thus,  is  fundamental  to  the 
deployment  of  the  technique  as  a  diagnostic. 

Under  Army  Research  Office  sponsorship.  Contract  DAAG29-C-83-0001 ,  UTRC 
has  Investigated  the  phenomenon  of  pressure-induced,  i.e.  collisional,  nar¬ 
rowing  in  nitric  oxide  and  water  vapor.  Well  controlled  laboratory  experi¬ 
ments  provided  data  for  comparison  with  predictions  from  model  calculations. 

NO  CARS  spectra  have  been  obtained  over  the  range  1-33  atmospheres,  300-800  K; 
H2O  CARS  spectra  from  1-30  atmospheres,  500-1000  K.  A  computer  code  was 
assembled  to  predict  these  CARS  spectra  and  comparisons  between  theory  and 
experiment  have  been  made. 

It  is  also  of  interest  to  Investigate  the  capability  of  CARS  for 
quantitative  measurements  of  the  density  and  temperature  of  chemical  species 
in  high  pressure  combustion  environments.  Quantitative  measurements  were 
studied  in  a  premixed  CO/air  flame  operated  at  pressures  up  to  15  atmospheres. 
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N2  and  CO  CARS  spectra  were  obtained  from  this  flame  as  well  as  Raman 
scattering  data  for  N2,  CO  and  C02«  The  experimental  data  were  compared 
to  chemical  equilibrium  code  predictions  and  predictions  from  the  UTRC  CARS 
code  for  the  flame  conditions.  During  the  course  of  this  contract,  a 
technical  advance  in  CARS  diagnostics  was  achieved  at  UTRC  through  the 
development  of  dual  broadband  CARS.  This  technique  overcomes  one  of  the 
previous  disadvantages  of  CARS,  namely  Interrogation  of  one  chemical  species 
at  a  time.  With  dual  broadband  CARS  a  multiplicity  of  species  can  be  sampled 
simultaneously,  greatly  expanding  the  utility  of  CARS  as  a  diagnostic. 

In  the  next  section  of  the  report,  the  spectroscopic  investigations  of 
H2O  and  NO  are  presented;  the  succeeding  section  describes  the  high  pressure 
flame  results  and  is  followed  by  a  section  on  dual  broadband  CARS.  Finally 
the  references  are  listed  and  appendices  are  attached  that  list  personnel 
involved  in  this  work,  publications/presentations  supported  under  this 
contract  and  a  reprint  of  a  recent  paper  on  dual  broadband  CARS. 


R85-956328-F 


HIGH  PRESSURE  CARS  SPECTROSCOPY 


Introduction 


A  coherent  anti-Stokes  Raman  vibrational  Q-branch  spectrum  Is  composed  of 
a  large  number  of  transitions  which  are  strongly  broadened  by  collisions  and 
the  molecule's  dipolar  intermolecular  potential.  At  high  pressures  the 
closely  spaced  lines  broaden  Into  one  another  and  the  spectrum  undergoes  a 
transition,  referred  to  as  colllslonal  narrowing.  In  which  spectral  features 
appear  to  coalesce.  Under  these  conditions,  the  usual  calculation  of  a  syn¬ 
thetic  spectrum,  carried  out  as  If  It  were  composed  of  Isolated  lines  Is  no 
longer  appropriate.  Beginning  In  1979  under  ARO  sponsorship,  UTRC  has  been 
studying  this  phenomenon  In  a  number  of  molecules  of  diagnostic  relevance  to 
the  Army.  Under  the  previous  contract,  N2  and  CO2  were  Investigated 
(Stufflebeam,  et  al . ,  1984;  Hall  and  Stufflebeam,  1984)  and  under  the  current 
effort  such  studies  have  been  extended  to  H2O  and  NO.  CO  was  studied  under 
an  Investigation  supported  by  the  Air  Force  Rocket  Propulsion  Laboratory 
(Stufflebeam,  et  al . ,  1983). 

In  the  earlier  contract.  Hall  et  al.  (1980)  developed  the  exact  CARS 
spectral  Intensity  expression  in  the  colllslonal  narrowing  regime  within  the 
Impact  approximation.  This  description  requires  a  knowledge  of  state-to-state 
relaxation  rates  and  Inversion  of  the  G-matrlx  at  every  point  In  frequency 
space.  The  technique  of  Gordon  and  McGinnis  (1968),  requiring  just  a  single 
matrix  dlagonallzatlon  and  Inversion  per  spectrum,  was  only  very  recently 
applied  by  Koszykowski  et  al.  (1985)  to  nitrogen  CARS  calculations.  Thus, 
this  approach  was  not  available  during  the  present  Investigations.  The  algo¬ 
rithm  represents  a  significant  computational  economy  and  will  be  Implemented 
Into  the  UTRC  CARS  codes  In  the  near  future. 

Because  of  the  large  number  of  transitions  possible  in  water  vapor, 
calculations  within  this  G-matrlx  framework  are  prohibitive  even  on  a  fast 
minicomputer,  and  therefore  approximate  formulations  have  been  sought.  Hall 
and  Greenhalgh  (1982)  extended  Gordon's  rotational  diffusion  concept  to 
describe  colllslonal  narrowing  In  N2  CARS  spectra  and  compared  the  model 
with  measurements  In  pure  water  vapor  at  pressures  up  to  16  atmospheres 
(Greenhalgh,  et  al.,  1984).  These  measurements  have  been  extended  In  this 
report  to  Include  nitrogen  broadening,  typical  of  actual  combustion 
situations.  The  CARS  spectrum  of  water  vapor  has  been  studied  in  flames  by 
Hall  et  al.  (1979),  at  lower  temperatures  in  a  heated  cell  with  a  number  of 
collision  partners  at  UTRC  (Hall  and  Shirley,  1983)  and  at  elevated  pressures 
at  AERE  Harwell  (Greenhalgh,  et  al.,  1984).  The  measurements  reported  here 
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were  undertaken  to  test  calculations  at  elevated  temperatures  typical  of 
practical  devices  with  air-fed  combustion.  In  such  devices  the  water  mole 
fraction  can  vary  from  a  fraction  of  a  percent  to  15-20  percent.  The 
measurements  undertaken  simulate  these  conditions  In  a  closed  cell  with  a 
nitrogen  collision  partner,  at  temperatures  from  500  to  1000  K.  These 
measurements  therefore  complement  the  Harwell  Investigations  and.  In  addition, 
the  spectral  resolution  In  the  reported  measurements  Is  higher,  posing  a  more 
critical  test  of  theoretical  calculations. 

The  CARS  spectrum  of  MO  has  been  studied  up  to  2  atmospheres  by  Laane  and 
Kiefer  (1980)  and  at  lower  pressure  by  Lempert,  et  al.  (1984).  This  report 
contains  measurements  of  the  CARS  spectrum  of  NO  at  pressures  up  to  33  atmos¬ 
pheres,  nearer  the  practical  combustion  environments  typical  of  propellants 
where  this  molecule  Is  an  Important  Intermediate  product  that  may  control 
energy  release  mechanisms. 

In  this  section  the  treatment  of  CARS  calculations  at  high  pressure  Is 
described,  followed  by  a  brief  description  of  pertinent  aspects  of  H2O  and 
NO  spectroscopy,  llnewldth  measurements  and  the  llnewldth  models  used  to  pro¬ 
vide  the  Input  for  computer  calculations.  The  experiments  are  then  discussed 
and  comparisons  are  made  with  theoretical  predictions.  Finally  the  conclu¬ 
sions  of  these  studies  and  directions  for  future  Investigations  are  discussed. 


Theoretical  Description  of  CARS  at  High  Pressure 


The  theoretical  basis  for  the  calculation  of  CARS  spectra  In  the  pressure 
narrowing  regime  Is  well  described  In  the  literature  (Hall,  et  al.,  1980;  Hall 
and  Greenhalgh,  1982;  Hall  and  Eckbreth,  1984;  Koszykowskl,  et  al.,  1985). 

The  purpose  of  this  section  Is  to  give  a  physical  picture  of  narrowing  and  to 
present  the  various  CARS  models  developed,  so  their  differences  can  be 
understood. 

The  frequencies  of  rotational  components  of  the  vibrational  Q-branch,  for 
which  there  Is  no  change  In  rotational  quantum  number  during  the  transition, 
differ  because  of  the  effect  of  rotation  on  the  vibrational  frequency,  l.e. 
vibrational-rotational  Interaction.  Therefore,  depending  upon  rotational 
state,  molecules  radiate  at  slightly  different  vibrational  frequencies, 
leading  to  a  simple  broadening  of  components  of  the  Q-branch.  When  the 
collision  frequency,  which  leads  to  Individual  line  broadening.  Is  small 
compared  to  the  vibration-rotation  frequency  shift,  the  Q-branch  Is  well 
resolved.  If  the  molecular  density  Is  Increased,  and  therefore  collisions 
occur  more  rapidly,  a  point  Is  reached  when  the  collision  frequency  is 
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comparable  to  the  frequency  shift  due  to  the  vibration-rotation  Interaction. 
The  Q  components  then  overlap  and  rotationally-inelastic  collisions  cause 
transitions  so  quickly  that  It  Is  possible  to  speak  only  of  the  aggregate  of 
molecules  In  an  average  J  state.  As  a  result  the  band  collapses  to  a  position 
determined  by  the  most  probable  rotational  state.  Note  that  broadening  can 
still  occur  through  vlbratlonally-lnelastic  and  dephaslng  collisions. 

The  Isolated  lines  approximation  is  valid  at  low  pressures.  The  CARS 
susceptibility  expression  then  is  just  the  result  of  simple  summation  of 
Individual  transitions: 


(3)  ^  27V 


(0) 


2(u}if  —  cjj  +  W2)  —  iT,-^ 


+  X 


Ew + ■)<;') 


where  N  is  the  molecular  number  density,  is  the  matrix  element 
of  the  polarizability,  is  the  density  matrix  element,  is  the 

Raman  transition  frequency,  ^If  is  the  Raman  linewidth  and  and  U2  are 

the  CARS  pump  and  probe  frequencies,  respectively.  X  is  the  nonresonant 
susceptibility.  Given  linewidths,  this  expression  is  straightforward  to 
evaluate. 


When  the  collision  frequency  is  high,  the  susceptibility  must  be 
evaluated  from  the  expression  derived  from  the  polarizability  and  impact 
approximations  (Hall  et  al.,  1980;  Gordon,  1966b). 

(w)  =  ^  ^  Or  a,  (2) 


where  the  elements  of  the  G-matrlx  are  given  by 

Grs  =  »(wi  -  -  Wr)^r«  +  [(rr/2)  —  lAr]^r»  +  (3) 

Here  the  y  are  the  components  of  the  matrix  elements  of  the  polarizability, 
is  the  frequency  shift  with  pressure,  6^^  is  the  Kronecker  delta  and 
the  are  the  of f-dlagnonal  linewidth  parameters,  describing  the  rate  of 
colllslonal  energy  transfer  between  states.  Calculation  of  CARS  spectra 
involves  Inversion  of  the  G-matrix  at  each  frequency,  u,  a  cumbersome  and  time 
consuming  chore  for  molecules  with  many  transitions,  like  water.  The  Gordon- 
McGinnis  technique  is  attractive  because  it  is  exact  and  only  one  matrix  diag- 
onalization  and  inversion  is  required  per  spectrum  however  it  was  unavailable 
during  this  study. 
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An  approximate  approach,  based  on  Gordon's  (1966a)  J-diffusion  model 
results  in  a  simpler  expression  to  evaluate; 


(3)  _  _2 


EAp;;V(2Au;j  -tTy) 

3 


l  +  *E(l-«^))r;P5;V(2Au;y-tTy) 


(4) 


where  cfi.  represents  the  fractional  contribution  of  pure  dephaslng  processes 
to  the  llnewidth  of  the  transition  and  Awj  is  the  detuning  frequency  factor 
seen  in  parentheses  in  the  denominator  of  the  first  expression  of  the  right 
hand  side  of  Eqn.  (1).  The  advantage  of  this  expression  is  that  it  requires 
little  more  time  to  evaluate  than  the  Isolated  lines  expression. 


For  moderate  pressures,  where  the  effects  of  overlap  are  small  G~^  can 
be  expanded  to  first  order  in  density  to  give  the  first  order  perturbation 
approximation: 


y(3)_2ArrY^  I  g.-y 

h  [^2AUr  -iVr  ^2Au3r-Yr^  2Au3s-Vs  ^  ^ 

Note  that  overlap  effects  appear  as  a  correction  to  the  isolated  lines 
approximation . 

In  calculating  the  CARS  spectrum  of  molecules  from  these  formulae,  two 
things  are  needed:  (1)  knowledge  of  the  spectroscopy  and  (2)  linewidths  (or 
equivalently  the  energy  transfer  rates).  The  spectroscopy  affects  energy 
levels,  transition  frequencies  and  polarizability  matrix  elements  (selection 
rules).  This  is  very  briefly  discussed  in  the  next  section. 


Water  Vapor 


Spectroscopy  of  Water  V'apor 

The  water  molecule  is  an  asymmetric  top.  It  has  three  fundamental 
vibrational  modes  corresponding  to  symmetric  stretching,  bending  and 
asymmetric  stretching  motions;  denoted  respectively  by  vibrational  quantum 
numbers  v^  V2  v^.  All  modes  are  Raman-active  but  the  vj  mode  at  a 
Raman  shift  of  3657  cm“^  is  by  far  the  strongest.  At  low  temperatures  only 
the  ground  state  is  significantly  populated,  and  Raman  transitions  belonging 
to  the  000  ->  100  band  are  most  important.  The  first  level  in  the  bending 
mode  has  an  energy  of  1597  cm“^,  so  the  vibrational  hot  band,  010  ■*  110, 
and  others  become  important  as  the  temperature  increases. 
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The  rotational  substructure  is  not  amenable  to  representation  by  an 
explicit  expression.  The  structure  Is  represented  by  the  rotational  quantum 
number  J  and  pseudo  quantum  number  t.  t  takes  on  values  In  the  range  -J<t<J. 
Rotational  energies  In  the  2J+1  sub-levels  Increase  as  t  Increases  from  -J  to 
+J. 


Modeling  the  transitions  present  In  the  CARS  spectrum  of  H2O  has  been 
described  In  many  papers  (Hall  and  Shirley,  1983;  Greenhalgh,  et  al.,  1984; 
Gordon,  1966b).  The  same  prescription  has  been  followed  here.  The  energy 
level  data  compiled  by  Bribes  et  al.  (1976)  and  supplied  by  R.  L.  St.  Peters 
has  been  used.  The  Raman  cross  section  for  water  Is  subject  to  uncertainty 
(Greenhalgh,  et  al.,  1984).  A  value  of  1.4  x  10~^®cm^/sr  has  been  used 
here.  Likewise  the  nonresonant  susceptibility  of  water  is  not  well  known,  a 
value  2.0  times  that  of  nitrogen  was  used  (Hall  and  Shirley,  1983). 

Water  Llnewldths 


Modeling  of  llnewldths  Is  crucial  In  the  formulation  of  CARS  spectroscopy 
of  H2O  at  high  pressures  (Greenhalgh,  et  al.,  1984).  In  the  pressure  broad¬ 
ened  regime  the  llnewldth  Is  governed  by  Inelastic  lifetime-limiting  colli¬ 
sions  and  by  elastic  collisions  which  dephase  the  vibrational  and  rotational 
motions.  Inelastic  collisions  In  which  a  change  occurs  In  vibrational  or 
rotational  energy  can  both  contribute.  Water  has  a  strong  dipole  moment  so 
rotatlonally  Inelastic  collisions  would  be  expected  to  be  Important.  Hall  and 
Eckbreth  (1984)  have  pointed  out  that  when  lifetime-limiting  rotatlonally 
Inelastic  collisions  predominate,  the  llnewldths  should  not  have  a  significant 
dependence  on  the  radiative  vibrational  and  rotational  selection  rules. 
Therefore  the  widths  for  microwave.  Infrared  and  vibrational  Raman  transitions 
are  expected  all  to  be  roughly  the  same.  This  Is  fortunate  because  the 
llnewldths  of  microwave  (Benedict  and  Kaplan,  1939;  1964)  and  Infrared 
transitions  have  been  the  subject  of  Investigations  (Mandln,  et  al.,  1982). 

Benedict  and  Kaplan  computed  llnewldths  of  water  microwave  transitions 
within  the  Anderson  theory  framework,  which  models  collisions  of  molecules  as 
following  straight  line  paths  under  the  Influence  of  an  attractive  potential. 
The  short  range  repulsion  Is  handled  with  a  'catch-all'  distance  of  closest 
approach  parameter.  The  transfer  of  rotational  energy  Is  controlled,  like 
many  energy  transfer  processes,  by  the  resonance  energy  concept,  which 
considers  the  energy  defect  In  the  collision  and  radiative  process.  The 
energy  defect  must  be  made  up  by  translational  energy,  the  availability  of 
which.  Is  of  course  determined  by  the  temperature. 

Benedict  and  Kaplan  (1964)  modeled  water-water  broadening  with  a  solely 
dlpole-dlpole  attractive  potential.  The  broadest  lines  they  calculated  were 
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transitions  involving  small  J  and  t.  Over  the  range  300<  T  <360  K,  the 
dependence  of  average  llnewldth  on  temperature  went  as  T~",  where  n  “  0.9. 

They  noted  that  extrapolation  of  the  widths  to  higher  temperatures  encountered 
In  flames  and  exhaust  was  not  warranted  for  individual  lines,  but  presumed 
that  the  average  n  for  entire  bands  will  approach  n  =  1.0,  expected  for 
resonant  dipole  collisions.  The  agreement  with  microwave  data  for  llnewldths 
at  a  particular  J  averaged  over  t  is  fairly  good. 

Benedict  and  Kaplan  (1959)  also  reported  calculations  of  nitrogen 
broadened  microwave  widths.  A  dlpole-quadrupole  attractive  potential  was 
used.  Calculated  widths  varied  from  0.2  to  0.06  cm“^  at  300  K,  with  the 
broader  widths  occurring  for  smaller  J  and  t.  Over  the  limited  temperature 
range  from  220  to  300  K  the  temperature  exponent  varied  widely  among  lines, 
with  n  =  0.8,  being  typical.  At  higher  temperatures,  however  no  such  simple 
variation  was  found  due  to  the  fact  that  at  higher  temperatures  all  lines  show 
a  degree  of  resonance. 

More  recently,  the  self  broadened  widths  of  Infrared  water  transitions 
have  been  reported  by  Mandln  et  al.  (1982).  They  calculated  widths  taking 
into  account  dipole-dipole,  dlpole-quadrupole,  quadrupole-dlpole  and 
quadrupole-quadrupole  Intermolecular  Interactions  within  the  theory  of 
Anderson-Tsao-Curnutte  (Tsao  and  Curnutte,  1962).  For  broadening  of  the  Vj 
transitions  these  and  the  Benedlct-Kaplan  calculations  are  comparable. 

The  greatest  Impediment  to  H2O  CARS  modeling  may  be  that  llnewldth  data 
at  high  temperature  is  lacking.  In  order  to  model  broadening  of  transitions 
at  high  temperatures,  the  exponential  gap  model  of  Polyani  and  Woodall  (1972) 
was  adopted.  In  this  concept,  explicit  selection  rules  governing  the  change 
of  rotational  quantum  number,  AJ,  are  suppressed,  but  the  rate  is  taken  to 
depend  exponentially  on  the  energy  defect.  This  approach  has  the  advantage  of 
simplicity  while  retaining  a  phenomenologically  correct  character.  The 
llnewldth  is  given  by: 

=  pao[T)Y^p^j\,txp{-c\t^Ejr,j>T'\lkT)  (6) 

J‘t> 


where  p  is  the  pressure,  a  parameter  corresponding  to  the  broadening 
coefficient  at  zero-energy  defect,  AEj^  j't'  energy  defect  in  the 

transition  and  c  is  a  fitting  parameter.  Tlie  parameter  Is  taken  -to  be  a 
function  of  temperature,  of  the  form  T“**.  The  n's  for  water  and  nitrogen 
broadening  are  taken  from  the  sensitivity  calculations  of  Benedict  and  Kaplan. 
The  other  parameters  have  been  found  by  fitting  llnewldths  In  the  exponential 
gap  model  to  the  measurements  of  Mandln  et  al.  (1982)  for  H2O-H2O,  and 
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Benedict  and  Kaplan's  (1959)  calculations  for  H2O-N2.  The  agreement  is 
shown  in  Fig.  1.  The  linewidths  shown  for  a  particular  J  have  been  averaged 
over  all  t  .  While  the  agreement  for  average  width  is  good  for  both  water  and 
nitrogen  broadened  widths,  if  individual  transitions  are  considered,  it  is 
seen  that  the  Benedict-Kaplan  calculations  tend  to  predict  a  broader  width  for 
small  1^1,  and  a  more  narrow  width  for  larger  1t|  .  This  could  be  an  Important 
discrepancy  because,  many  of  the  peaks  observed  in  H^O  CARS  spectra  corre¬ 
spond  to  the  overlap  of  the  more  highly  populated  1'^]  J  transitions.  The 
parameters  assumed  for  the  calculations  shown  are  summarized  in  Table  1. 


Table  1. 

Exponential  Gap  Llnewldth  Model  Parameters 


H2O-H2O  H2O-N2 


1.025  0.20 

n  1.0  0.8 

c  0.08  0.25 


Nitric  Oxide 


Spectroscopy  of  NO 

The  electronic  structure  of  the  NO  molecule,  more  fully  discussed  in 
Lempert  et  al.  (1984),  is  complicated  by  ground  state  splitting  that  results 
in  the  overlap  of  two  Q-branch  series  in  the  CARS  spectrum.  The  two  series 
are  separated  by  approximately  0.2  cm~^  and  the  two  series  are  not  resolved 
by  the  equipment  used  in  the  current  measurements  (0.5  cm“^).  Collisional 
coupling  between  the  rotational  levels  of  the  two  states  has  initially  been 
ignored  in  the  computer  model  but  may  be  incorporated  in  a  later  revision,  if 
needed. 

NO  Linewidths 


The  linewidths  used  in  the  NO  CARS  code  are  the  recent  measurements  of 
Lempert,  et  al.  (1984).  For  NO,  the  number  of  rotational  transitions  is 
tractable  and  use  of  the  full  matrix  inversion,  Eqns.  (2)  and  (3),  is 
preferable  to  the  Gordon  model  of  rotational  diffusion  because  it  was  shown  to 
be  more  accurate  in  N2  and  CO  (Stuf flebeam,  et  al.,  1983;  Stufflebeam,  et 


Linewidth  (FWHM),  cm  yatm 
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al.,  1984;  Koszykowski,  et  al.,  1985).  The  pressure  broadening  of  the  lines 
is  modeled  by  a  form  of  the  Inverse  power  law,  previously  used  successfully  in 
N2.  The  Inverse  power  law  is  employed  to  model  the  off  diagonal  matrix 
elements  >^.5! 


nrs  =  -Kopi°MErer‘^  (7) 


These  are  related  to  the  theoretical  llnewidths  through  detailed  balance 
and  conservation  of  probability. 


(8) 


The  theoretical  llnewidths  are  then  least-squares  fit  to  the  experimental 
ones  to  determine  the  fitting  parameters  K^,  a. 


a  a 

^0'  a-i 


=  0 


(9) 


This  procedure  produces  the  fit  (lines)  to  experimental  widths  shown  in 
Fig.  2.  The  Inverse  power  law  model  does  not  fit  the  measured  llnewidths 
very  accurately.  The  discrepancy  is  largest  at  low  J  for  the  n  »  3/2  state. 
It  is  expected  from  previous  work  in  CO  at  this  laboratory  (Stuff lebeam,  et 
al.,  1983)  that  another  rate  model  originally  proposed  by  Dexhelmer,  et  al. 
(1982) 


Ira  —  -Ko(2r+  (10) 


or  the  exponential  gap  model,  Eqns.  6  and  7,  will  fit  the  experimental  line- 
widths  more  closely.  Evaluation  of  these  models  in  the  NO  CARS  code  will  be 
the  next  step  in  evaluation  of  che  code. 
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H2O  Experimental  Measurements 


Experimental  Apparatus 

The  CARS  apparatus  used  In  these  experiments  has  been  described 
previously  (Hall  and  Shirley,  1983).  The  laser  beams  focus  through  the  water 
vapor  In  a  colllnear  CARS  arrangement.  The  Stokes  wave  Is  provided  by  a 
broadband  dye  laser  operated  with  either  LD-690  dissolved  In  methyl  alcohol  or 
DCM  dissolved  In  dimethyl  sulfoxide  (DMSO),  (Eckbreth  and  Anderson,  1985). 

The  DCM  mixture  was  centered  better  with  respect  to  the  Stokes  band  and  had 
higher  output  power. 


The  high  pressure  cell  used  In  these  experiments  Is  shown  In  Fig.  3.  It 
Is  made  up  of  heated  cell  within  an  outer,  pressure-containing  cell.  The 
outer  cell  has  been  hydrostatically  tested  to  2200  psl,  but  Is  not  meant  to 
withstand  high  temperatures.  The  water  vapor  Is  mostly  contained  within  the 
Inner  cell  which  Is  sealed  off,  except  for  the  connection  to  the  water 
reservoir.  The  Inner  cell,  therefore  Is  pressurized  by  nitrogen  through  the 
reservoir.  A  small  diameter  vent  tube  prevents  water  from  being  forced  out  of 
the  reservoir  as  the  pressure  rises  In  the  Inner  cell  due  to  heating.  The 
temperatures  of  the  reservoir  and  Inner  cell  are  controlled  separately.  In 
operation,  the  temperature  of  the  Inner  cell  Is  brought  up  to  the  desired 
operating  point  and  the  temperature  In  the  reservoir  Is  slowly  raised,  varying 
the  water  concentration  within  the  Inner  cell.  To  prevent  condensation  on  the 
outer  cell  windows,  the  area  surrounding  the  reservoir  Is  filled  with 
molecular  sieves. 

From  a  knowledge  of  the  reservoir  temperature,  the  concentration  of  water 
vapor  In  the  Inner  cell  can  be  calculated  in  theory  from  water  p-V-T  data.  It 
appears,  however,  because  the  cell  Is  not  in  equilibrium,  that  this  approach 
Is  not  reliable.  The  concentration  of  water  vapor  in  the  Inner  cell  Is  not 
well  known,  therefore.  There  are  no  Instruments  capable  of  measuring  water 
vapor  concentrations,  or  relative  humidity  In  the  extreme  environment  of  the 
cell. 


It  was  also  found  that  the  laser  flux  through  the  Inner  windows  had  to  be 
limited  to  prevent  optical  damage,  because  the  windows  were  so  close  to  the 
focal  point.  Sapphire  windows,  expected  to  have  superior  chemical  Inertness 
with  respect  to  the  high  temperature,  high  pressure  water  vapor,  could  not 
withstand  as  high  a  laser  flux  as  quartz  (probably  due  to  lower  quality 
surface  finish  on  the  sapphire).  To  obtain  the  highest  quality  spectra, 
quartz  windows  were  used.  Degradation  was  not  a  serious  problem  for  the 
limited  time  which  the  windows  were  exposed. 
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H2O  CARS  Experiments 

Measurements  were  undertaken  In  a  matrix  of  pressure  and  temperature 
combinations.  The  Intended  measurements  ranged  over  the  limits  of  normal 
operating  pressures  In  combustion  devices,  and  temperatures  attainable  in  the 
heated  cell.  Measurements  were  made  at  500,  750  and  1000  K  at  pressures  of  5, 
10,  and  20  atmospheres.  The  initial  pressure  In  the  cell  was  set  to  these 
values  and  was  permitted  to  rise  as  the  sealed  off  cell  was  heated.  Once  the 
cell  test  temperature  was  reached,  spectra  were  recorded  as  the  reservoir 
temperature  was  raised,  increasing  the  concentration  of  water  in  the  Inner 
cell. 


Figures  4  and  5  show  CARS  spectra  of  water  vapor  measured  at  500  K  for 
several  different  total  pressures.  Figure  4  Is  spectra  of  low  water  vapor 
concentration  and  Fig.  5  Is  of  high  concentration  spectra.  The  experimental 
spectral  resolution  Is  about  2  cm~^  for  these  data.  The  spectra  are  shown 
displaced  vertically  by  0.2  units  for  clarity.  The  modulation  due  to  the 
nonresonant  background  Is  clearly  apparent  in  Fig.  4.  Features  are  more 
sharply  defined  In  the  low  pressure  spectrum  as  expected.  This  is  true  In  the 
high  concentration  spectra  as  well,  (Fig.  5).  The  estimated  concentration  in 
Fig.  5  is  15-20%,  while  in  Fig.  4  the  concentration  is  probably  about  5%. 

CARS  spectra  of  high  water  vapor  concentrations,  measured  at  750  and 
1000  K,  are  shown  in  Figs.  6  and  7,  for  successively  higher  pressures.  At  the 
lowest  pressure  and  750  K  the  spectral  atures  are  quite  distinct,  but  lose 
contrast  at  the  higher  pressures.  The  same  is  true  at  1000  K.  Note  that 
whereas  some  of  the  features  of  the  band  closest  to  the  bandhead  are  somewhat 
resolved  at  500  K  at  the  highest  pressure,  at  the  higher  temperatures  the  same 
band  is  totally  unresolved  at  high  pressure.  At  the  highest  pressures  the 
cell  heater  was  unable  to  attain  1000  K.  Consequently  the  actual  measured 
temperature  at  29  atmospheres  was  931  K  and  at  12  atmospheres  it  was  985  K. 

Results  and  Discussion 


Model  calculations  have  been  prepared,  guided  by  the  results  of  the 
experiments.  Calculations  with  full  matrix  Inversion  have  not  been  done, 
because  of  the  prohibitive  time  required  for  calculations  with  the  full  set  of 
energy  levels  (831)  on  a  VAX/750  with  a  floating  point  accelerator.  On  the 
other  hand,  It  was  found  that  the  rotational  diffusion  and  perturbation  models 
gave  similar  results  v.hen  similar  linewldth  parameters  were  used. 

The  Isolated  lines  and  rotational  diffusion  calculations  at  500  K 
showed  no  dlscernable  difference  at  6  atoiospheres.  At  13  atmospheres  the 
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features  remained  more  distinct  in  the  rotational  diffusion  model,  for  the 
same  assumed  concentration.  The  diffusion  model  agreed  with  experiments,  in 
showing  a  rapid  fall  in  intensity  away  from  the  bandhead  on  the  low  energy 
side,  in  contrast  to  isolated  lines  model  calculations.  At  23  atmospheres  the 
low  energy  tall  was  about  50%  higher  for  Isolated  lines.  This  difference 
could  cause  a  considerable  error  in  concentration  measurements. 

The  experimental  trends.  Illustrated  in  Figs.  4-7,  have  been  examined 
with  the  rotational  diffusion  model.  Calculations  at  temperatures  and  pres¬ 
sures  comparable  to  the  experiments  are  shown  in  Figs.  8-10.  The  general 
trends  are  certainly  predicted  by  the  calculations,  which  used  the  linewidth 
parameters  in  Table  1,  with  no  contribution  for  dephaslng  collisions.  Clearly 
I  however,  at  lower  pressures  spectral  features  are  more  distinct  in  the  experi¬ 

ments  than  the  calculations  show.  Reducing  the  parameter  to  one-half  of 
the  Mandin  et  al.  (1982)  and  Benedlct-Kaplan  (1964)  values  resulted  in  sharper 
features  at  500  K,  but  at  1000  K  such  rates  show  the  spectral  features  as 
remaining  distinct,  in  disagreement  with  the  measurements.  Adding  a  dephasing 
contribution  may  affect  this,  but  there  has  been  no  corroborating  evidence  for 
*  such  a  scheme . 

Figure  11  shows  a  direct  comparison  of  model  calculations  and  an  experi¬ 
ment  in  the  middle  of  the  experimental  range.  The  experimental  data  are  shown 
as  points  and  the  full  line  gives  the  calculation.  The  assumed  water  concen- 
I  tratlon  has  been  adjusted  to  fit  the  band  head  and  the  tail  region.  The  ex¬ 

perimental  spectral  features  are  more  distinct  than  the  calculation.  The 
spectral  bandwidth  used  in  the  calculations  has  also  been  adjusted  to  fit  the 
experiment,  but  it  is  found  that  unreasonable  values  are  required  with  the 
linewidth  parameter  set  in  Table  1  to  obtain  a  fit  to  the  experiments. 

If  one  examines  the  transitions  which  cause  the  appearance  of  the  spec¬ 
trum  an  interesting  observation  is  revealed.  First  consider  the  grouping 
of  features  in  the  spectrum.  The  first  four  features  arise  from  the  funda¬ 
mental  band  000  -*•  100.  The  contribution  of  these  levels  is  falling  rapidly, 

when  at  about  3640  cm“^,  the  first  hot  band,  010  -*■  110,  arises.  The  next 

I 

major  group  is  the  second  hot  band,  020  -*■  120.  Looking  closely  at  transi¬ 
tions  within  these  groups,  they  appear  to  correspond  in  many  cases  to  the 
overlap  of  transitions  arising  with  |t|  >  J.  These  are  the  transitions  that 
have  larger  widths  in  the  exponential  gap  model  than  appear  reasonable  from 
Mandin  et  al.  and  Benedict-Kaplan's  calculations.  Smaller  linewidths  would 
I  sharpen  these  features,  to  give  the  higher  contrast  present  in  the  experi¬ 

mental  spectra. 
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Conclusions  and  Recommendations 

It  appears  at  present  that  the  collision  narrowing  models  used  in  CARS 
calculations  are  not  the  prime  limitation  to  water  vapor  calculations.  Rather 
it  seems  to  be  more  the  case  that  the  deficiency  is  in  the  detailed  behavior 
of  the  llnewldths.  Before  progress  can  be  made,  an  independent  means  of  de¬ 
termining  the  concentration  is  necessary.  Unfortunately  no  known  instrumen¬ 
tation  exists  which  is  suitable  for  application  in  the  relevant  parameter 
range.  The  best  approach  may  be  to  measure  absorption  on  a  well  resolved 
infrared  transition  with  a  tunable  diode  laser.  A  fruitful  avenue  of  research 
may  be  that  of  incorporating  llnewidth  calculations  into  the  CARS  code. 

Measurements  of  widths,  with  high  resolution  inverse  Raman,  for  example, 
would  be  very  useful.  Calculations  show  Important  differences  among 
transitions,  including  alternating  widths.  It  is  expected  however  that  it 
would  not  be  possible  to  resolve  all  important  lines,  because  of  the  combined 
width  and  proximity  of  lines  with  large  |t|. 

Future  improvements  include  use  of  the  Kataoka  et  al.  (1982)  and  Teets 
(1984)  convolution,  which  is  more  appropriate  for  measurements  with  current 
laser  sources.  It  is  is  expected  that  these  calculations  would  improve  the 
fit  of  spectra  near  the  bandhead.  Recently  Koszykowskl  et  al.  (1985)  reported 
the  use  of  the  Gordon  and  McGinnis  (1968)  technique  with  CARS  calculations  to 
eliminate  repeated  matrix  inversions  needed  to  calculate  a  complete  spectrum. 
With  this  simplification  a  matrix  needs  to  be  inverted  only  once. 

It  is  clear  that  further  work  is  necessary  before  CARS  can  be  used  for 
concentration  determinations  of  water  in  high  pressure  combustion  situations. 


NO  Experimental  Measurements 


NO  CARS  spectra  have  been  obtained  over  the  range  1-33  atmospheres, 
300-800  K.  Specific  combinations  of  these  parameters  were  selected  for  the 
acquisition  of  CARS  spectra  as  indicated  In  Fig.  12.  A  computer  code  was 
assembled  to  predict  the  NO  CARS  spectra  and  comparisons  between  theory  and 
experiment  have  been  made 

Experimental  Apparatus 

To  provide  well  controlled  conditions  of  temperature  and  pressure  for  the 
acquisition  of  NO  CARS  spectra,  the  high  pressure  vessel  described  in 
Stuff lebeam,  et  al.  (1984)  was  utilized.  It  is  a  static,  internally  heated 


Temperature  (K) 


R85-956328-F 


R85-956328-F 


cell  that  Is  rated  for  temperatures  to  1750  K  and  pressures  to  5000  psig.  It 
Is  equipped  with  a  pressure  gauge,  type  S  thermocouples,  and  an  optical 
pyrometer  to  measure  the  pressure  and  temperature  of  the  sample.  Only  the 
central  15  cm  Is  heated  and  sapphire  rods  extend  Into  the  heated  region  to 
reduce  thermally-induced  density  gradient  effects  in  the  gas.  The  limited 
optical  access  provided  by  the  test  cell  dictates  either  colllnear  alignment 
of  the  pump  and  Stokes  lasers  used  to  generate  CARS  or  a  tightly  packed 
BOXCARS  arrangement.  A  folded  BOXCARS  phase  matching  geometry  was  employed 
for  the  NO  measurements  with  an  axial  resolution  of  approximately  4  mm  to 
avoid  any  problems  from  axial  temperature  gradients  in  the  cell.  The  UTRC 
High  Pressure  CARS  Facility  provided  the  hardware  for  these  measurements  and 
is  fully  described  in  the  experimental  apparatus  section  of  the  high  pressure 
flame  measurements  of  this  report. 

NO  Experiments 


A  problem  was  encountered  during  the  NO  investigations  that  prohibited 
acquisition  of  NO  spectra  at  all  potential  conditions  of  interest  (see  Fig. 
12).  NO  is  highly  susceptible  to  catalytic  reactions  for  temperatures  below 
1000  K  and  this  presented  a  problem  for  acquisition  of  CARS  spectra  from  the 
high  pressure  cell,  ^s  mentioned  in  Stufflebeam,  et  al.  (1984)  the  cell 
contains  many  refractory  oxide  materials  such  as  alumina  and  zirconla  and  its 
heater  is  pure  platinum.  All  these  surfaces  are  highly  reactive  with  NO  as 
investigated  many  years  ago  by  Wise  and  Freeh  (1955),  Fraser  and  Daniels 
(1958),  Yuan,  et  al.  (1959),  among  others.  The  catalytic  reaction  converts 
the  NO  to  NO2  which  is  highly  absorbing  for  visible  wavelengths.  Calcula¬ 
tions  indicate  only  0,5%  NO2  at  1  atmosphere  pressure  will  produce  a  50% 
reduction  in  the  expected  CARS  signal.  The  reduction  accounts  for  absorption 
of  the  pump  laser  in  addition  to  the  absorption  of  the  CARS  frequency.  Tests 
of  this  effect  were  performed  by  loading  the  heated  cell  with  100%  NO  and 
monitoring  the  absorption  of  the  beam  that  was  propagated  through  the 
gas.  At  1  atmosphere  pressure,  the  NO  would  allow  transmission  of  the  pump 
beam  and  remain  optically  clear  for  more  than  one  and  one-half  hours  up  to  a 
temperature  of  800  K.  At  300  K  the  cell  remained  optically  cleai;  to  the 
highest  available  pressure  of  33  atmospheres  for  extended  periods  also.  The 
CARS  data  acquisition  mode  utilized  for  the  NO  spectra  employs  broadband 
generation  and  scanned,  narrowband  detection.  Since  the  scanning  CARS  tech¬ 
nique  requires  up  to  40  minutes  to  obtain  the  CARS  spectrum,  and  time  is 
required  for  thermal  equilibrium  after  filling  the  cell  with  NO,  it  was  judged 
that  any  observable  NO2  generation  in  less  than  one  hour  would  jeopardize 
the  experiment.  Our  absorption  tests  indicated  this  was  possible  at  1  atmos¬ 
phere  from  300-800  K  and  at  300  K  from  1-33  atmospheres  (Fig.  12).  For 
pressures  above  1  atmosphere  the  usable  temperature  range  was  so  limited,  up 
to  about  430  K,  as  to  preclude  useful  information  to  be  gained  for  this 
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investigation.  The  above  limits  were  also  confirmed  through  an  analysis  using 
the  CHEMKIN  chemical  kinetics  code  to  predict  NO2  production  using  activa¬ 
tion  energies  and  frequency  factors  for  the  heterogeneous  reaction  published 
by  Wise  and  Freeh  (1955)  with  additional  reactions  published  by  McCullough,  et 
al.  (1977).  There  was  insufficient  time  to  design  and  fabricate  a  new  high 
pressure  cell  with  materials  that  would  not  react  with  NO,  so  only  those 
spectra  indicated  in  Fig.  12  were  acquired. 

CARS  spectra  of  NO  were  acquired  at  295  K  over  the  pressure  range  1-33 
atmospheres  and  the  bandwldths  (FWHM)  measured  from  the  experimental  data. 
These  data  are  presented  in  Fig.  13,  where  the  different  symbols  represent 
data  acquired  on  different  days.  As  with  other  molecules  studied  at  UTRC, 
there  is  little  change  in  bandwidth  as  the  pressure  is  varied  and  it  is 
obvious  the  colllsional  narrowing  process  is  operative  over  this  pressure 
range . 


The  experimental  spectra  acquired  at  1  atmosphere,  616  and  801  K  are 
presented  in  Figs.  14  and  15.  The  abrupt  discontinuity  in  these  spectra  in 
the  region  of  20,620  cm~^  is  artificially  produced  during  the  scan  to  check 
the  zero  signal  level.  It  is  observed  that  even  the  616  K  spectrum  exhibits  a 
hot  band  and  that  both  spectra  show  significant  modulation  of  the  nonresonant 
background  in  the  vicinity  of  the  fundamental  bandhead.  The  nonresonant 
background  susceptibility  of  NO  is  3  times  that  of  N2  so  modulation  is 
expected  at  these  temperatures  even  for  high  concentration  of  the  active 
molecule.  Computer  predictions  are  in  progress  for  these  experimental 
conditions. 

Results  and  Discussion 


Shown  in  Fig.  16  is  the  experimental  CARS  spectrum  of  NO  at  295  K,  1 
atmosphere  (dots)  together  with  the  prediction  of  the  NO  CARS  code  for  the 
same  conditions.  The  computer  prediction  used  the  measured  NO  llnewidths  of 
Lempert  et  al  (1984)  and  th-  inverse  power  law  model  to  determine  the  off- 
diagonal  llnewldth  elements  o'  the  G  matrix  (Stuf flebeam,  et  al,  1984).  Other 
models  such  as  isolated  lines  or  the  Gordon  model  of  rotational  relaxation 
also  predict  this  case  well,  but  the  Inverse  power  law  reproduces  the  spectra 
most  closely  over  the  whole  pressure  range.  Figure  17  is  the  same  as  Fig.  16 
except  the  pressure  is  20  atmospheres.  The  fit  is  again  very  good  and  for 
this  case  the  other  models  of  llnewldth  do  not  compare  as  favorably  with  the 
experiment. 
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Recommendations 

The  basis  for  the  NO  predictive  CARS  code  is  the  laser  bandwidth  con¬ 
volution  proposed  by  Yuratich  (1979)  to  model  the  effect  of  finite  bandwidth 
pump  lasers  on  the  predicted  CARS  spectrum.  In  recent  years  this  approach  has 
been  questioned  by  researchers  in  Japan  (Kataoka,  et  al.,  1982)  and  in  the 
U.S.  (Teets,  1984).  They  proposed  a  new  convolution  that  has  been  shown, 
through  rigorous  comparisons,  to  give  better  agreement  with  experimental  data. 
It  is  known  as  the  Kataoka-Teets  (K-T)  convolution  and  is  discussed  more  fully 
in  connection  with  high  pressure  flame  measurements  later  in  this  report.  The 
improvements  are  most  pronounced  for  spectra  of  resonant  species  in  low  con¬ 
centration  and  when  the  ratio  of  pump  laser  width  to  Raman  linewidth,  Aw/F,  is 
large.  It  is  not  expected  to  have  a  large  effect  for  the  current  studies  be¬ 
cause  they  deal  with  high  concentration  of  the  active  molecule  and  the  param¬ 
eter  Acj/r  gets  smaller  at  high  pressure.  It  would  appear  necessary,  however, 
to  convert  the  current  NO  CARS  code  to  one  based  on  the  K-T  convolution  for 
predictions  of  NO  spectra  at  low  concentration  and  low  pressure  and  this 
effort  is  currently  in  progress. 

Future  directions  of  this  work  Involve  Implementation  of  the  Dexhelmer 
power  law  (Eqn.  10)  and  exponential  gap  (Eqn.  6)  models  of  rotational  relaxa¬ 
tion  in  the  NO  CARS  code  as  well  as  to  fit  the  1  atmosphere  high  temperature 
data.  It  is  expected  from  previous  results  (Stuff lebeam,  et  al.,  1983)  that 
the  Dexhelmer  power  law  will  predict  the  NO  CARS  spectrum  more  closely  over 
the  entire  range  of  temperature  and  pressure  studied.  It  will  also  be 
interesting  to  compare  calculations  using  the  exponential  gap  model  to  the 
results,  as  Koszykowskl,  et  al.  (1985)  achieved  good  comparisons  with  this 
model  for  N2  CARS  spectra  up  to  10  atmospheres.  To  decide  on  which  of  these 
models  is  best,  comparisons  are  needed  with  high  pressure,  hot  spectra.  These 
studies  will  be  difficult  unless  special  cells  are  built  to  prevent  the  NO 
NO2  conversion  problem  encountered  here. 
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HIGH  PRESSURE  FLAME  MEASUREMENTS 


Introduction 


As  part  of  the  current  contract,  UTRC  extended  the  previous  studies  of 
high  pressure  CARS  spectra  of  various  molecules  to  the  regime  of  high  pressure 
flames.  For  these  studies  a  CO/alr  premixed  flame  was  chosen  because  of  the 
relatively  simple  chemistry  Involved.  Furthermore  there  are  no  hydrocarbons 
or  soot  formed  which  can  complicate  interpretation  of  temperature  and  concen¬ 
tration  measurements  from  the  CARS  spectra  (Dobbs  et  al.,  1985).  A  premixed 
flame  was  chosen  to  avoid  the  complication  of  having  to  account  for  diffusion 
effects  In  the  flame.  The  products  of  the  flame  chemistry  and  their  relative 
concentrations  form  a  useful  base  for  CARS  measurements.  Nitrogen  Is  the 
major  species  and  Is  used  for  thermometry.  CO  Is  present  In  low  enough 
concentration  to  demonstrate  concentration  measurements  by  spectral  shape. 

Also  as  will  be  shown,  both  spectra  can  be  obtained  simultaneously  with  the 
apparatus  employed  for  this  work.  This  turns  out  to  be  advantageous  for 
temperature  and  concentration  measurements  as  will  become  clear  In  later 
sections  of  this  report. 

The  next  section  describes  the  high  pressure  flame  facility  developed  for 
this  contract.  Then  the  CARS  and  Raman-scatterlng  apparatus  are  described. 
Finally,  the  results  of  the  CARS  and  ivaman  studies  of  the  high  pressure  flame 
are  presented  and  discussed;  directions  for  future  work  are  also  Indicated. 


Experimental  Apparatus 


High  Pressure  Flame  Facility 

The  premixed  flame  was  chosen  to  operate  In  non-sootlng  fashion  and  with 
a  fairly  simple  chemical  system.  Consequently,  the  simple  fuel  CO,  was 
selected.  The  choice  of  CO  as  the  fuel  gave  a  very  simple  product  gas  compo¬ 
sition  (CO,  CO2,  N2,  and  possibly  O2)  to  expedite  atom  balance  closure 
for  qtiantltatlve  CARS  evaluation.  This  avoided  complicating  the  CARS  measure¬ 
ments  with  the  Inclusion  of  H2,  H2O,  and  other  hydrogenlc  species.  The 
oxidation  of  CO  to  CO2  Involves  an  Intermediate  reaction  with  an  OH  radical. 

It  was  necessary  to  Include  a  trace  amount  of  a  gas  that  would  promote  the 
Ignition  and  reaction  of  CO  with  air.  Consequently,  ultra  pure  CO  Is  used  as 
the  fuel  and  a  trace  amount  of  C2H^  Is  Introduced  Into  the  fuel  stream  as 
a  source  of  the  OH  Intermediate.  This  allowed  the  flame  to  be  well  controlled 
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over  the  entire  range  of  operating  pressure.  It  also  allowed  some  control 
over  the  flame  speed. 

A  high  pressure  test  cell,  designed  and  constructed  earlier  was  employed 
to  house  the  burner.  The  vessel  is  a  heavy  wall  cylinder  with  top  and  bottom 
flanges  sealed  with  0-rings.  Optical  access  is  provided  by  four  large  aper¬ 
ture  (2.5  cm)  quartz  windows  at  90  degrees,  thus  permitting  CARS,  LIF  or  Raman 
measurements  to  be  made.  Provision  has  also  been  made  for  thermocouples, 
pressure  transducers,  gas  supply,  ignition,  and  exhaust.  The  flameholder, 
designed  for  flexibility  of  application,  is  replaceable  and  can  be  changed  by 
removing  the  bottom  flange.  The  inside  surface  of  the  quartz  windows  are 
located  about  6  cm  from  the  centerline  of  the  test  cell  to  prevent  damaging 
the  windows  by  the  incident  high  energy  laser  beams. 

A  premixed  flame  burner  was  designed,  constructed,  and  employed  to 
investigate  quantitative  CARS  measurements  in  a  high  pressure  CO/air  flame. 

The  burner  design  uses  a  1.3  cm  diameter  fuel  tube  capped  with  a  brass  tip 
that  has  37  drilled  holes  (1.0  mm  diameter)  in  a  closely-packed  hexagonal 
pattern.  Figure  18  illustrates  the  burner  design  and  shows  the  details  of  the 
brass  burner  tip.  The  flame  was  operated  with  the  cavity  beneath  the  burner 
tip  filled  with  No.  12  grade,  HP  copper  shot,  0.15  cm  in  diameter.  The  feed 
tube  is  surrounded  by  a  3.8  cm  o.d.  (0.09  cm  wall)  stainless  tube  that 
provides  the  secondary  shroud  flow,  generally  nitrogen.  There  was  no 
observable  difference  in  the  premixed  flame  whether  or  not  the  shroud  flow  was 
used.  Consequently,  the  shroud  flow  was  not  used  in  the  CARS  experiments. 

The  flame  is  ignited  at  atmospheric  pressure  with  a  spark  from  a  removable 
wire  centered  over  the  burner  tip  and  driven  by  a  commercial  Tesla  coll.  The 
Ignitor  wire  is  positioned  over  the  center  of  the  premixed  burner  so  that  a 
spark  of  about  4-5  mm  in  length  is  obtained. 

Provision  was  made  to  have  five  gases  available  to  the  burner;  carbon 
dioxide,  nitrogen,  air,  carbon  monoxide,  and  ethylene.  Standard  gas  cylinders 
are  connected  to  Hastings  mass  flowmeters  and  then  to  needle  valves  which 
provide  the  flow  control.  The  flowmeters  were  calibrated  with  a  wet  test 
meter.  The  two  fuel  lines,  CO  and  C2H^,  and  the  air  line  are  connected  to 

a  mixing  chamber,  thence  to  the  fuel  tube  of  the  burner.  The  mixing  chamber 

cavity  is  filled  with  No.  9  lead  shot  which  distributes  and  mixes  the  flow. 

The  pressure  in  the  high  pressure  test  cell  is  controlled  by  regulating 
the  exhaust  flow  out  of  the  cell  through  a  needle  valve.  The  flow  rates  are 
adjusted  to  maintain  a  stable  flame  at  the  desired  stoichiometry  and  as  the 
cell  is  pressurized,  the  exhaust  needle  valve  is  adjusted  to  maintain  the  test 

cell  pressure.  Once  stable  conditions  have  been  reached,  the  cell  operates 

for  hours  without  any  apparent  change  In  burner  characteristics  or  operating 
pressure. 
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FIG  18 
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The  brass  tip  with  the  small  holes  stabilizes  the  flame  front  about  1  mm 
above  the  surface.  The  premixed  flame  burner  was  operated  successfully  at 
pressures  up  to  about  15  atmospheres,  providing  a  flat  premixed  flame  about  1 
cm  In  diameter  at  atmospheric  pressure  and  a  slightly  reduced  diameter  at 
elevated  pressure  as  shown  in  Fig.  19.  It  was  necessary  to  Increase  the  flow 
rate  of  fuel  and  air  to  the  burner  at  elevated  pressure  to  prevent  quenching 
of  the  flame  on  the  burner  surface.  It  was  also  possible  to  blow  the  flame 
off  the  burner  If  the  feed  rate  was  too  high. 

Experiments  were  performed  to  determine  the  practical  operating  envelope 
of  the  burner  at  various  flow  rates,  air/CO  ratios,  and  pressure  levels 
(Hedman,  et  al.,  1986).  A  set  of  "design"  gas  flow  rates  was  selected  for 
CARS  measurements  at  each  operating  pressure.  The  design  air/CO  ratio 
selected  Is  about  1.6,  chosen  to  reside  well  within  the  stable  operating 
envelope.  Based  on  adiabatic  chemical  equilibrium  calculations  (Gordon  and 
McBride,  1976),  these  operating  conditions  give  the  product  gas  distribution 
and  adiabatic  flame  temperature  listed  in  Table  2. 


Table  2 

Gordon-McBride  Predictions  of  CO/alr  Premixed  Flame 


Pressure 

%  CO 

%  CO^ 

7, 

AFT 

1  atm 

18.5 

26.1 

53.4 

2295  K 

5 

16.9 

27.2 

54.2 

2334 

10 

16.3 

28.0 

54.4 

2352 

CARS  Apparatus 


The  experiments  were  carried  out  in  the  UTRC  High  Pressure  CARS  Facility 
shown  schematically  In  Fig.  20.  The  output  of  a  Quanta-Ray  NdrYAG  laser  is 
frequency  doubled  to  generate  a  horizontally  polarized,  "primary"  pump  beam  at 
532  nm  (u^).  The  10  Hz  laser  has  an  average  output  of  2.2  Watts  at  532  nm 
with  a  pulse  duration  of  7  nanoseconds.  A  beamsplitter  separates  30Z  from 
1^2  CO  pump  the  broadband  Stokes  dye  laser  oscillator  and  amplifier.  The 
primary  pump  is  subsequently  passed  through  a  polarizer  then  divided  Into  two 
pump  beams  for  the  folded  BOXCARS  phase-matching  geometry  (Shirley,  et  al., 
1980)  by  a  beamsplltter/mlrror  combination.  The  polarization  of  each  beam  is 
Independently  controlled  by  half-wave  plates.  An  optical  delay  line  Is 
Inserted  in  one  of  the  legs  to  establish  the  degree  of  correlation 
desired  for  the  two  pump  beams.  An  optical  flat  is  inserted  after  the  Stokes 
beam  amplifier  to  displace  ^2  plane  of  the  two  primary  pumps. 
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FIG  20 
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Galilean  telescopes  are  provided  to  control  beam  waists  and  the  focal  zone 
locations  of  u 2  and  the  wj's.  CARS,  at  Is  generated  In  the  focal 

volume,  and  all  three  frequencies  (wj,  602,  103)  are  recollimated  after 
exiting  the  cell.  The  CARS  beam  is  separated  from  the  pump  beams  by  a  spatial 
filter  then  passed  through  a  broadband  polarization  analyzer,  for  the 
background  suppression  experiments  (Eckbreth  and  Hall,  1981),  before  the 
signal  is  incident  on  the  slit  of  a  0.5  m  spectrograph.  A  He-Ne  laser  is 
shown  whose  output  Is  coincident  with  the  (^2  and  used  for  optical 

alignment.  A  beamsplitter  In  the  1U3  leg  provides  a  reference  signal  to  a 
photomultiplier  which  is  used  to  monitor  the  spectrally-integrated  CARS  signal 
strength.  A  broadband  diode  array  detector  Is  mounted  on  the  spectrograph. 

It  Is  a  1000  channel,  intensified  array  that  allows  acquisition  of  very  broad 
CARS  spectra  ('-600  cm”^).  Use  of  the  vidlcon  greatly  expedites  signal 
collection  and  reduces  measurement  times  considerably.  The  data  (spectra) 
from  the  diode  array  are  sent  to  the  computer  for  spectral  processing  and 
archival  storage. 

Raman  Apparatus 

To  complement  CARS  measurements  in  the  high  pressure,  premixed,  CO/alr 
flame,  a  Raman  scattering  diagnostic  capability  was  designed  and  assembled  for 
the  UTRC  High  Pressure  CARS  Facility.  The  Raman  diagnostic  system  shown  In 
Fig.  21,  consists  of  two  photomultiplier  tube  detectors,  monitored  through 
gated  integrators  for  single  laser  pulse  temporal  response.  Their  output  is 
recorded  by  A/D  converters  In  a  laboratory  computer.  One  tube  detects  Stokes- 
shifted  Raman  scattering  from  the  532  nm  pump  laser  and  the  other  tube  detects 
anti-Stokes  shifted  photons  from  the  same  laser  pulse.  The  laser  polarization 
was  rotated  to  vertical  for  the  Raman  measurements.  Several  narrowband  inter¬ 
ference  filters  were  designed  and  obtained  so  that  various  flame  species  (CO, 
CO2,  N2)  could  be  monitored.  The  Raman  scattered  radiation  was  imaged  on 
an  aperture  to  control  the  size  of  the  observed  experimental  volume.  Addi¬ 
tional  shields  and  absorbing  glass  filters  were  employed  to  reduce  interfer¬ 
ence  from  scattered  radiation  at  532  nm.  The  system  was  calibrated  for  each 
filter  by  Raman  scattering  from  a  gas  at  known  conditions  (l.e.  N2  at  NTP). 
Data  from  the  flame  then  yielded  the  number  density  (Stokes  shifted  scatter¬ 
ing)  of  the  observed  gas  species  in  the  combustion  environment.  Similarly, 
the  anti-Stokes  radiation  was  calibrated  by  a  tungsten  lamp  at  known  tempera¬ 
ture  and  the  experimental  data  was  recorded.  The  ratio  of  anti-Stokes  to 
Stokes  scattering  gives  an  estimate  of  the  gas  temperature  (Eckbreth,  at  al., 
1979). 


Photographs  of  the  experimental  apparatus  are  shown  in  Figs.  22  and  23. 
The  high  pressure  combustion  vessel  and  Raman  apparatus  are  evident  in  Fig. 
22;  the  pump  laser  and  beam  manipulation  optics  used  to  generate  the  folded 
BOXCARS  phase  matching  geometry  are  shown  in  Fig.  23. 
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Experimental  Measurements 


Raman  Measurements 


Using  the  previously  described  Raman  apparatus,  measurements  of  product 
gas  relative  concentration  were  taken.  The  flame  was  stabilized  at  each 
operating  pressure  and  the  output  of  the  gated  integrator  that  monitors  the 
Stokes  Raman  PMT  was  sampled  by  an  A/D  converter.  Five  hundred,  single-pulse 
datapolnts  were  obtained  for  each  gas  species  of  interest.  Different  gases 
were  measured  by  changing  the  narrow  band  ( ~ 40  cm“^  FWHM)  interference 
filter  in  front  of  the  PMT.  Each  filter  was  designed  to  minimize  its 
temperature  sensitivity  vis-a  vis  the  number  of  resonant  species  photons 
collected  (Eckbreth,  at  al.,  1979).  In  this  way  the  collected  signal  is 
proportional  to  the  number  density  of  the  scattering  species  and  independent 
of  the  temperature.  Each  filter  combination  was  calibrated  at  room  tempera¬ 
ture  with  the  cell  filled  with  100%  of  the  respective  gas  that  was  being 
measured.  The  calibrations  then  are  used  to  relate  measured  signals  to  gas 
number  density  and  relative  concentrations.  The  gas  temperatures  can  be 
inferred  through  Raman  measurements  described  in  Eckbreth  et  al.  (1979).  The 
basic  technique  is  to  measure  the  ratio  of  anti-Stokes  to  Stokes  scattering 
for  the  same  molecule,  in  this  case  nitrogen.  The  Raman  apparatus  was 
designed  to  obtain  this  measurement  and  the  data  were  acquired  for  each 
operating  pressure  of  the  flame.  At  this  time  the  temperature  data  have  not 
been  fully  analyzed;  however,  the  analysis  is  ongoing  and  results  are  expected 
in  the  near  future.  The  Raman  measurements  of  relative  species  concentrations 
in  the  premixed  flame  are  presented  in  Fig.  24.  The  symbols  Indicate  experi¬ 
mental  measurements  and  the  lines  are  the  predictions  from  the  Gordon-McBride 
computer  code.  The  agreement  for  N2  and  CO2  is  within  7%  and  quite  good. 

The  comparisons  for  CO  are  less  favorable  (20%)  and  Indicate  a  need  for 
further  analysis.  Since  the  experimental  data  are  uniformly  higher  than  the 
computer  prediction,  a  recallbratlon  of  the  CO  Stokes  detector  may  be 
necessary. 

CARS  Measurements 


The  UTRC  high  pressure  flame  facility  provided  the  medium  to  obtain  CARS 
spectra  from  the  high  pressure  premixed  flame.  The  CARS  spatial  resolution 
was  measured  as  a  .16  mm  dla.  x  3.5  mm  long  cylinder  which  is  adequate  to 
probe  the  1  cm  dia.  high  pressure  flame.  The  pump  beams  were  crossed  in  the 
folded  BOXCARS  phase  matching  geometry  approximately  2  mm  above  the  burner 
surface.  This  is  above  the  reaction  zone  so  that  only  flame  products  were 
measured.  Examples  of  CARS  spectra  from  various  operating  pressures  of  the 
flame  are  shown  in  Figs.  25-28.  All  of  these  data  are  from  the  CO-N2 
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Raman-resonant  region  and  effects  of  pressure  are  most  noticeable  as  the 
narrowing  and  sharpening  of  the  fundamental  bandhead  of  the  CARS  signatures. 
Previous  work  supported  by  ARO  (Contract  DAAG29-79-C-0008)  has  shown  that 
colllslonal  narrowing  Is  the  dominant  mechanism  operational  in  this  pressure 
range.  The  operating  pressure  limit  of  the  current  premixed  flame  Is 
obviously  adequate  to  produce  the  colllslonal  narrowing  phenomenon  in  the  CARS 
signatures  so  that  higher  operating  pressures  were  not  sought. 


Results  and  Discussion 


The  approach  to  CARS  measurements  In  the  high  pressure  flame  Involves 
analyzing  the  CARS  spectral  distribution  (i.e.  intensity  variation  with 
frequency)  and  Is  fully  described  In  Hall  and  Eckbreth  (1984).  A  CARS  spec¬ 
trum  such  as  Fig.  25  Is  composed  of  resonant  and  nonresonant  contributions. 

The  resonant  parts  (I.e.  CO  and  N2)  are  functions  of  the  resonant  molecule 
number  density,  temperature,  resonant  molecular  constants,  and  pressure, 
through  the  rovlbratlonal  llnewldths  of  the  molecule.  The  nonresonant  contri¬ 
bution  Is  dispersionless,  a  function  of  total  gas  density  and  Implicitly  tem¬ 
perature,  through  the  Ideal  gas  law.  Assuming  the  pressure  Is  Independently 
measured  and  accurate  values  of  molecular  constants  are  available,  data  such 
as  Fig.  25  contain  Information  on  the  concentrations  of  CO,  N2,  the  total 
gas  density,  and  the  gas  temperature.  The  difficulty  in  extracting  the 
information  lies  in  separating  the  various  functional  dependencies  to  yield 
estimates  of  the  experimental  parameters.  In  general,  any  point  of  the 
Intensity  distribution  of  Fig.  25  is  dependent  to  some  degree  on  all  the  above 
parameters.  In  mathematical  terms  the  data  Is  represented  by  the  squared 
modulus  of  the  third  order,  nonlinear  susceptibility  (Eqn.  1), 

Ixl^  =  X'^  +  2x'x'‘’’ +  X”’’*  +  X"^  (11) 


where  x'  and  x"  are  the  real  and  Imaginary  parts  of  the  resonant  suscepti¬ 
bility  and  x*^*^  la  the  nonresonant  part.  Examination  of  Eqn.  (11)  can  guide 
data  analysis.  For  instance,  it  Is  clear  that  even  when  there  Is  no  resonant 
response  In  the  sampled  frequency  interval  i.e.  x',  x"  “  0»  there  is  still  a 
CARS  signal  represented  by  x°^  which  is  proportional  to  total  gas  density 
and,  by  invoking  the  ideal  gas  law,  temperature.  If  a  molecule  has  a 
resonance  in  the  frequency  window  but  is  in  low  concentration  such  that 
(xtesonantj2  ^^nrj2^  ^  modulation  of  the  nonresonant  spectrum  is 
predicted  through  the  2x  ”’^x'  term.  For  this  low  density  region,  the 
spectral  shape  is  density  and  temperature  dependent  and  potentially  useful  for 
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measurements  of  these  parameters.  CARS  Is  a  unique  spectroscopy  in  this 
regard,  concentration  can  be  measured  from  spectral  shape  in  certain  density 
ranges.  The  technique  is  limited  to  relatively  low  concentration,  say  less 
than  30%,  and  yields  only  a  relative  concentration  to  the  total  density 
represented  by  the  nonresonant  part.  For  higher  concentration  of  a  resonant 
species  the  spectral  shape  concentration  sensitivity  is  lost  and  concentration 
measurements  must  rely  on  an  absolute  measurement.  For  molecules  in  very  high 
concentration,  the  background  (nonresonant)  signal  is  lost  and  the  pure 
resonant  response  is  most  evident.  This  response  is  very  sensitive  to 
frequency  and  temperature,  and  of  course  the  integrated  resonant  response  is 
proportional  to  concentration  of  the  resonant  molecule  squared,  although  the 
frequency  distribution  of  the  resonant  part  is  not  as  concentration 
sensitive  at  this  high  level.  This  region  of  high  density  is  generally 
utilized  for  thermometry  because  the  CARS  intensity  distribution  with 
frequency  is  largely  just  a  function  of  temperature  if  the  resonant  gas 
concentration  is  high. 

The  foregoing  ideas  have  been  incorporated  in  computer  models  that 
predict  CARS  spectra  given  the  total  gas  temperature  and  density  and  the 
concentration  of  the  resonant  species.  By  comparing  an  experimental  spectrum 
against  a  library  of  computed  predictions,  a  CARS  signature  can  be  'least 
squares'  fit  for  temperature  and  species  concentration;  the  reliability  of  the 
estimates  depends  on  the  concentration  range  of  the  data.  Generally  tempera¬ 
ture  is  most  reliable  from  the  CARS  spectrum  of  a  species  in  high  concentra¬ 
tion  and  this  is  why  N2  thermometry  is  usually  employed.  Estimates  of  the 
concentration  from  spectral  fits  are  most  reliable  for  concentrations  in  the 
range  1-30%. 

UTRC  has  developed  CARS  predictive  codes  and  fitting  routines  based  on 
equations  published  by  Yuratlch  (1979).  The  equations  account  for  the 
convolution  of  pump  laser  llnewidth  and  resonant  Raman  linewidth  to  produce 
the  observed  CARS  frequency  distribution.  In  general  the  solutions  are  very 
complex  but  have  been  solved  for  some  cases  of  practical  Interest.  The 
solution  employed  here  assumes  Gaussian  laser  statistics.  Measurements  were 
conducted  using  the  Quanta  Ray  DCR-IA  laser  of  the  UTRC  high  pressure  CARS 
facility  to  investigate  the  mode  statistics  of  the  pump  laser  in  a  manner 
similar  to  the  experiments  of  Rahn  et  al.  (1984).  Using  the  folded  BOXCARS 
phase  matching  geometry,  an  optical  delay  was  introduced  in  one  leg  of  the 
primary  pump.  With  the  laser  operating  at  a  linewidth  of  ~  0.8  cm"^,  the 
correlation  length  of  the  optical  pulse  Is  ~  1.2  cm.  The  spectrally 
integrated  CARS  signal  generated  in  both  a  nonresonant  gas  and  resonant  N2 
in  air  at  295  K  were  monitored  on  a  photomultiplier/gated  integrator  detector. 
The  two  signals  and  their  ratio,  as  a  function  of  the  optical  delay,  are 
presented  in  Fig.  29.  Following  the  analysis  of  Rahn  et  al.  (1983),  the  data 
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describe  a  laser  with  Gaussian  mode  statistics.  Of  significance  for  our  CARS 
experiments  and  In  contrast  to  the  Sandla  experiments  Is  that  the  ratio  of 
resonant  to  nonresonant  signals  Is  Independent  of  the  pump  beam  correlation. 

It  Is  an  advantage  to  operate  at  full  correlation  (0  delay)  since  the  signal 
strengths  are  highest  and  the  signal  to  noise  ratio  Is  maximized. 

Using  the  UTRC  fitting  code,  which  utilizes  the  Yuratlch  convolution  the 
data  of  Fig.  25  were  analyzed  for  temperature  and  concentration  by  spectral 
shapes.  A  comparison  of  experimental  data  (dots)  and  computer  fit  (line)  Is 
presented  In  Fig.  30.  The  match  between  theory  and  experiment  Is  excellent 
for  a  theoretical  temperature  of  1922  K  and  concentration  of  38%  N2.  The 
concentration  measurement  Is  below  that  predicted  by  the  Gordon-McBrlde 
equilibrium  calculation  and  the  Raman  scattering  measurements  and  may  be  a 
result  of  the  Insensitivity  to  concentration  of  spectral  shape  measurements  In 
this  concentration  region.  The  predicted  temperature  Is  reasonable  since  the 
adiabatic  flame  temperature  Is  2295  K  and  we  expect  radiation  losses  and 
conduction  losses  to  the  metal  burner.  Using  the  same  code,  the  CO  portion  of 
Fig.  25  can  also  be  fit  for  temperature  and  resonant  species  concentration  and 
those  results  are  presented  In  Fig.  31.  The  fit  again  appears  to  be  excellent 
and  this  time  the  estimated  concentration  (17.2%)  Is  closer  to  the  equilibrium 
calculation  of  18.5%.  However  the  predicted  temperature  Is  not  the  same  as 
that  predicted  for  the  N2-resonant  segment  of  the  same  spectrum  and  moreover. 
Is  very  near  the  adiabatic  flame  temperature  which  is  physically  unrealistic. 
Resolution  of  this  dilemma  resides  In  the  laser  convolution  used  to  predict 
CARS  spectra  for  the  fitting  routine. 

During  the  period  of  this  contract  several  researchers  (Kataoka,  et.  al., 
1982;  Teets,  1984)  have  questioned  the  validity  of  the  Yuratlch  laser  convolu¬ 
tion  calculations.  Kataoka  and  later  Teets  presented  alternate  solutions  of 
the  convolution  equations  that  they  showed  to  be  more  accurate  for  certain 
experimental  conditions.  Farrow  and  Rahn,  (1985)  also  presented  data  to 
verify  the  new,  (K-T),  convolution.  When  the  pump  laser  linewldth  Is  very 
narrow  with  respect  to  the  Raman  transition  linewldth  of  Interest,  the 
Yuratlch  solution  Is  valid.  For  the  case  where  the  pump  laser  linewldth  Is 
larger  than  the  Raman  linewldth,  additional  terms  have  to  be  Included  In  the 
equation.  The  Quanta  Ray  laser  has  a  linewldth  of  ~  0.3  cm”^  at  532  mm  when 
used  with  the  optional  Intra  cavity  etalon  and  a  linewldth  of  ~  0.8  cm~^ 
when  the  etalon  is  not  used.  Either  condition  Is  larger  than  typical  Raman 
llnewldths  In  atmospheric  flames  ('-0.03  cm~^).  Greg  Dobbs  of  our  group  has 
programmed  a  CARS  predictive  code  based  on  the  K-T  convolution.  A  comparison 
of  this  computer  prediction  with  the  data  of  Fig.  25  is  shown  in  Fig.  32.  A 
rigorous,  least  squares  fit  was  not  possible  for  this  case  as  the  K-T  predic¬ 
tive  code  Is  undergoing  evaluation  and  has  not  yet  been  Incorporated  In  the 
CARS  fitting  routine.  The  results  of  Teets'  work,  confirmed  by  Farrow  and 
Rahn,  Indicate  the  major  Impact  of  the  K-T  convolution  is  for  spectra  of  low 
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concentration  gases,  and  that  temperature  predictions  based  on  the  K-T  convo¬ 
lution  for  majority  species  will  not  be  very  different  from  that  predicted  by 
the  Yuratich  convolution  for  the  same  majority  species.  Starting  from  the 
assumption  then,  that  the  predicted  N2  temperature  from  Fig.  30  is  nearly 
correct,  a  series  of  visual  fits  were  prepared  for  the  data  that  represents 
the  resonant  nitrogen  part  of  the  experimental  spectrum.  This  analysis  showed 
(Fig.  32)  that  the  concentration  predicted  by  the  chemical  equilibrium  code 
and  as  measured  by  the  Raman  scattering  data  fit  the  experimental  data 
closely.  Also,  variations  in  predicted  temperatures  were  evaluated  and  the 
1922  K  estimate,  based  on  the  Yuratich  result  was  best.  The  estimated 
temperature  and  concentration  were  1922  ±  25  K  and  53  ±5%,  based  on  this 
visual  fitting.  Figure  33  presents  a  similar  visual  fit  of  the  resonant  CO 
portion  of  the  same  experimental  data  (Fig.  25).  The  starting  assumption  for 
this  visual  fit  was  that  the  CO  temperature  should  be  the  same  as  the  N2 
temperature.  This  assumption  was  borne  out  in  the  fitting  procedure,  with  the 
additional  result  that  the  predicted  concentration  is  in  agreement  with  the 
Gordon-McBrlde  and  Raman  scattering  data. 

It  is  satisfying  to  achieve  the  foregoing  consistent  set  of  data  and 
computer  predictions.  It  appears  clear  from  this  work  that  the  K-T  convolu¬ 
tion  needs  to  be  Incorporated  in  the  CARS  spectral  fitting  routines.  A 
theoretical  aspect  still  needs  to  be  addressed  however;  the  K-T  convolution 
has  not  been  rigorously  proven  applicable  to  an  experimental  sitxjation  with 
correlated  pump  lasers.  In  the  work  at  Sandia  (Farrow  and  Rahn,  1985),  where 
the  lasers  apparently  exhibited  non-Gaussian  statistics,  decorrelation  of  the 
pump  beams  was  necessary  to  obtain  spectral  agreement  with  the  K-T  theory. 

With  the  apparent  Gaussian  statistics  here,  good  agreement  was  obtained  with 
correlated  pump  beams.  Further  investigation  is  necessary  in  this  area. 

Measurements  from  the  CARS  spectra  at  high  pressure  are  more  cumbersome 
than  the  1  atmosphere  case.  The  primary  reason  for  this  is  that  the  UTRC 
high  pressure  CARS  code  is  predictive  at  this  time  and  not  incorporated  into  a 
computer  fitting  routine.  Thus  visual  overlays  are  required  to  obtain 
estimates  of  temperature  and  concentration.  Based  on  the  previous  results  of 
the  1  atmosphere  case,  the  K-T  convolution  will  be  required  at  high  pressure. 
Although  the  ratio  of  laser  linewidth  Co  Raman  llnewldth  is  smaller  at  high 
pressure,  Che  ratio  should  still  be  in  the  region  where  K-T  is  required  i.e. 
Auj/r>  1.  The  K-T  convolution  has  not  yet  been  implemented  for  the  high 
pressure  case  at  UTRC  so  concentration  measurements  from  spectral  shapes  of 
high  pressure  CARS  spectra  such  as  Figs.  26-28  could  not  be  performed.  The 
resonant  nitrogen  portion  of  the  spectrum  should  still  yield  reliable 
temperature  measurements  through  analysis  with  the  Yuratich  laser  convolution, 
however.  It  was  shown  in  previous  work  at  UTRC,  (Stufflebeam,  et  al.,  1983; 
198A)  that  high  pressure  CARS  spectra  are  most  accurately  predicted  by  a  full 
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matrix  Inversion  model  of  the  resonant  signature.  Simpler,  faster  models  such 
as  Isolated  lines,  perturbation  expansion  or  the  Gordon  model  of  rotational 
relaxation  were  inadequate  for  the  pressures  encountered  in  these  experiments. 
The  need  to  revert  to  matrix  inversion  slows  the  analysis  Immensely.  Twenty- 
four  hours  of  CPU  time  on  a  VAX  750  are  required  to  produce  one  CARS  spectrum 
from  the  high  pressure  predictive  code  using  matrix  inversion  with  the  Inverse 
power  law  model  used  to  predict  the  off  diagonal  matrix  elements.  The  algo¬ 
rithm  of  Gordon  and  McGinnis,  (1968),  not  yet  Incorporated  in  the  UTRC  codes, 
should  help  this  situation.  The  Gordon  model  is  applicable  to  high  pressure 
and  much  faster  (^3  min.  of  CPU  time)  although  less  accurate.  It  was  used  to 
narrow  the  temperature  range  of  analysis  and  then  a  full  matrix  spectrum  was 
generated  for  confirmation  with  the  experimental  data.  The  result  for  the  10 
atmosphere  case  (Fig.  27)  is  shown  in  Fig.  34.  The  temperature  predicted  by 
the  visual  fitting  is  2150  K  and  seems  physically  resonable  based  on  the 
adiabatic  flame  temperature  of  2352  K.  The  'best  fit'  concentration  was  80% 
which  seems  clearly  in  error  and  reinforces  the  need  for  K-T  convolutions  at 
high  pressures.  Based  on  the  extensive  computer  time  required  for  these  fits 
and  the  suspicious  concentration  estimates  and  overpredictions  of  temperature 
derived  from  Yuratlch  convolutions  of  low  concentration  spectra,  an  attempt 
was  not  made  to  analyze  the  CO  portion  of  the  experimental  data  at  10 
atmospheres  (Fig.  27). 

Additional  techniques  can  be  used  to  measure  resonant  species  concentra¬ 
tion  from  CARS  data.  These  are  appropriate  when  the  density  is  high  enough  to 
invalidate  spectral  shape  fitting.  The  techniques  Involve  absolute  Intensity 
measurements  and  must  be  calibrated  to  a  known  set  of  conditions  of  density 
and  temperature.  They  require  monitoring  the  resonant  species  signal  intensi¬ 
ty,  Integrated  over  all  transition  frequencies  of  the  molecule,  and  have  to  be 
normalized  to  the  pump  laser  Intensities  I^l2« 

As  mentioned  previously  the  CARS  spectrum  generally  contains  contribu¬ 
tions  from  a  resonant  species  and  the  dispersionless,  nonresonant  background. 
When  the  resonant  species  signal  strength  is  much  greater  than  the  nonresonant 
contribution,  an  absolute  Intensity  measurement  merely  requires  splitting  off 
a  portion  of  the  total  CARS  signal  to  a  photomultiplier  and  calibrating  the 
signal  to  a  known  standard  gas,  in  effect  the  nonresonant  signal  is  swamped  by 
the  resonant  signal  and  neglected.  Pump  laser  normalization  could  be  achieved 
by  directly  measuring  the  intensities  on  a  photodiode  but  this  ignores  space 
and  time  effects  of  the  Intensity  in  the  experimental  volume  and  can  be  dan¬ 
gerous  for  all  but  the  most  benign  combustion  environments.  It  is  preferable 
to  generate  a  nonresonant  signal  in  an  exterior  cell  of  known  gas  concentra¬ 
tion  and  temperature.  The  resultant  calibration  signal  from  the  nonresonant 
cell  would  also  be  proportional  to  I^l2  could  be  used  for 
normalization. 
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The  more  general  approach,  and  one  necessary  for  the  high  pressure  CO/alr 
premixed  flame,  addresses  the  case  when  the  nonresonant  contribution  to  the 
CARS  spectrum  cannot  be  neglected.  In  this  Instance  the  resonant  and 
nonresonant  photons  need  to  be  separated  and  individually  measured.  The 
technique  is  called  polarization  suppression  and  is  fully  described  in  Hall 
and  Eckbreth  (1984).  The  resonant  and  nonresonant  signals  have  different 
responses  to  polarization  directions  of  the  laser  fields.  By  adjusting  the 
polarization  angles  of  the  pump  beams,  the  resonant  and  nonresonant  signals, 
which  are  generated  at  different  polarization  angles,  can  be  separated.  The 
purely  resonant  signal  is  then  monitored,  normalized  and  calibrated  to 
complete  the  quantitative  measurement.  The  normalization  and  calibration  can 
be  handled  in  a  'series'  or  'parallel'  nonresonant  cell  as  described 
previously  but  the  method  is  sensitive  to  alignment  errors  introduced  by 
refractive  index  differences  between  the  'experimental'  (flame)  path  and  the 
path  through  the  calibration  cell,  normally  held  at  NTP.  A  more  attractive 
normalization  is  provided  by  the  'suppressed'  nonresonant  signal  itself 
(Farrow,  et  al.,  1982;  1985).  By  monitoring  the  purely  nonresonant 
contribution  I^l2  normalization  is  achieved  and  an  'In-sltu'  calibration 
is  provided  because  the  nonresonant  signal  is  proportional  to  total  gas 
density  at  exactly  the  same  location  as  the  resonant  signal.  Scaling  of  the 
dispersionless,  nonresonant  signal  to  conditions  of  known  density  and 
temperature  are  accomplished  by  filling  the  reaction  vessel  with  a  known  gas 
at  a  known  temperature.  Being  suitably  calibrated,  the  nonresonant  signal  is 
then  the  calibration  standard  for  the  resonant  signal.  Providing  that  the 
nonresonant  susceptibility  for  the  experimental  gas  can  be  calculated  (usually 
a  valid  assumption),  and  that  the  temperature  is  measured  from  the  resonant 
species  signature,  the  total  gas  density  of  the  nonresonant  signal  provides 
the  calibration  for  the  resonant  response. 

The  polarization  suppression  techniques  are  not  widely  used  with 
broadband  approaches  because  they  invoke  at  least  a  factor  of  16  loss  of  the 
resonant  signal  (Hall  and  Eckbreth,  1984).  However,  for  high  pressure  regimes 
this  reduction  may  be  tolerable.  The  polarization  suppression  technique  was 
partially  Implemented  in  the  high  pressure  premixed  flame.  Initial  experi¬ 
ments  were  conducted  to  demonstrate  its  applicability  to  the  high  pressure 
regime.  For  this  experiment  half-wave  plates  rotated  the  plane  of  polariza¬ 
tion  of  each  wj  pump  beam  to  60°  relative  to  the  'P'  (0°)  polarization 
of  the  U2  Stokes  laser.  This  combination  of  angles  maximizes  the  separation 
of  the  resonant  and  nonresonant  parts  of  the  spectrum.  A  portion  of  the  pure 
resonant  contribution  is  viewed  through  an  analyzer  rotated  to  120° 
(perpendicular  to  the  nonresonant  signal).  The  result  for  the  10  atmosphere 
flame  is  shown  in  Fig.  35.  The  data  of  Fig.  35  exhibited  no  pressure  or 
temperature  Induced  depolarization  (Farrow,  R.  L.,  et  al.,  1984;  Ikawa  and 


R85-956328-F 


Whalley,  1984)  from  the  quartz  windows  of  the  combustion  vessel.  The  spectrum 
In  Fig.  35  demonstrates  very  good  background  suppression,  as  Indicated  by  the 
lack  of  modulation  at  the  bandheads  of  the  two  resonant  species  and  the 
zero-signal  baseline  (compare  with  Fig.  27). 

By  rotating  the  analyzer  to  150°  (perpendicular  to  the  resonant 
contribution)  a  portion  of  the  nonresonant  signal  Is  obtained  to  provide  an 
'In  situ'  normalization  for  laser  power  fluctuations  and  a  calibration  for 
concentration  measurements  of  the  resonant  species.  Further  work  will 
Incorporate  separate  photomultlpller/gated  integrator  detector  channels  for 
the  purely  resonant  and  nonresonant  signals  and  calibration  experiments  using 
a  gas  at  known  concentration  and  temperature.  Suitable  analyses  of  these 
signals  will  produce  data  for  the  gas  temperature  and  the  total  gas  and 
resonant  species  concentration. 


Recommendations 


These  research  Investigations  at  UTRC  are  continuing  and  future 
directions  Include  completion  of  the  partially  Implemented  N2  concentration 
measurements  and  Implementation  of  the  polarization  suppression  technique  for 
concentration  measurements  of  CO2  as  well.  The  spectroscopy  of  CO2  has 
been  studied  under  previous  Army  funding  (Contract  DAAG29-79-C-0008)  and  the 
results  of  the  CO2  Investigations  appeared  In  Hall  and  Stufflebeam,  (1984). 
Although  CO2  is  expected  to  be  In  low  concentration  in  the  high  pressure 
flame  the  very  small  spacing  of  rotational  transitions  In  the  CO2  CARS 
spectrum  results  In  significant  line  overlap  In  each  vibrational  band,  even  at 
1  atmosphere.  The  strong  overlap  Is  evidenced  by  a  colllslonally  narrowed 
spectrum  with  little  or  no  destructive  Interference  with  the  nonresonant 
background  (Hall  and  Eckbreth,  1984).  This  condition  Is  not  amenable  to 
concentration  measurements  from  spectral  shapes  and  so  the  more  complex 
technique  of  polarization  suppression  is  warranted. 

These  Initial  results  are  quite  encouraging  vis-a-vis  the  ability  of  CARS 
to  make  quantitative  measurements  In  high  pressure  combustion  environments. 

The  directions  are  clear.  The  K-T  laser  convolution  needs  to  be  incorporated 
In  the  predictive  computer  codes  to  Improve  the  accuracy  of  temperature  and 
concentration  predictions.  Additional  rotational  relaxation  models  (l.e. 
exponential  gap)  need  to  be  Incorporated  for  comparative  predictions  of 
colllslonal  narrowing  and  inclusion  of  the  Gordon-McGlnnls  algorithm  should 
greatly  speed-up  computer  response.  These  enhancements  of  the  diagnostic 
ability  of  CARS  should  contribute  to  programs  of  relevance  to  the  Army. 
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DUAL  BROADBAND  CARS 


During  this  contract  period,  a  novel  approach  to  multiple  species 
measurements  using  vibrational  CARS  In  which  only  one  additional  laser  beam 
needs  to  be  Introduced  was  developed  at  UTRC  (Eckbreth  and  Anderson,  1985, 
Appendix  C  of  this  report).  One  of  the  limitations  of  CARS,  as  normally 
implemented,  is  the  inability  to  interrogate  more  than  one  constituent  at  a 
time.  Except  for  a  few  fortuitous  Instances  of  closely  spaced  Raman 
resonances  (i.e.  N2  and  CO  as  evidenced  in  the  previous  section  of  this 
report),  a  separate  Stokes  dye  laser  Is  required  for  each  constituent  to  be 
measured.  Using  two  separate  Stokes  lasers,  N2  and  O2  have  been  simulta¬ 
neously  measured  (Jarrett,  et  al.  1985).  In  reality,  to  employ  more  than  two 
Stokes  sources,  becomes  quite  cumbersome  and  complex. 

The  dual  broadband  technique  employs  two  Independent  broadband  Stokes 
sources  In  combination  with  a  pump  laser  in  the  following  manner.  Besides  the 
two  separate  two-color  wave-mixing  processes  between  the  pump  and  Stokes 
lasers,  spectrally-resolved  CARS  Is  produced  in  a  three-color  process  from 
species  whose  Raman  resonances  correspond  to  the  frequency  differences  of  the 
two  broadband  sources.  Particularly  serendipitous  in  hydrocarbon-fueled 
combustion  is  the  fact  that  when  the  Stokes  sources  are  positioned  to  coincide 
with  CO2  and  H2O,  the  major  products  of  combustion,  the  frequency  differ¬ 
ence  between  the  two  coincides  with  the  Raman  resonances  of  heavy  dlatomlcs 
such  as  N2,  CO,  NO  which  can  then  be  observed  via  the  three-color  process. 
Yuratlch  (1979)  pointed  out  that  one  could  perform  CARS  with  high  spectral 
resolution  using  two  broadband  sources  If  the  third  laser  source  was  narrow- 
band.  Since  these  Raman  frequencies  are  well-defined,  spectral  smearing  does 
not  occur  even  with  broadband  excitation.  The  narrowband  laser  then  scatters 
off  the  coherent  Raman  excitation  In  the  medium  resulting  In  CARS  at  a  spec¬ 
tral  resolution  characteristic  of  the  narrowband  source. 

One  broadband  source  covers  the  1200-1600  cm”^  Stokes  spectral  region 
as  illustrated  In  Fig.  36.  This  source,  wavemixing  with  a  frequency-doubled 
neodymlum:YAG  laser  beam  at  532  nm,  produces  CARS  In  the  normal  way  from 
CO2,  O2  and  potentially  C-'’  bond  stretches  should  such  constituents  be 
present.  A  second  broadband  dye  laser  covers  the  3000-3700  cm“^  region  and 
produces  CARS  in  the  normal  way  from  H2O  and  C-H  stretch  Raman  modes.  The 
dye  to  use  is  DCM  which  lases  over  a  very  wide  spectral  region  around  640  nm. 
The  two  broadband  dye  lasers  in  turn  excite  all  Raman  resonances  with  shifts 
in  the  range  from  about  1400  cm“^  to  2500  cm”^  and  produce  CARS,  after 
wavemixing  with  532  nm,  from  such  constituents  as  N2,  N2O,  CO,  HCN,  NO  and 
CbC  stretches.  Despite  the  use  of  two  broadband  sources  in  the  three  color 
wavemixing,  spectral  resolution  is  determined  by  the  pump  laser  linewidth 
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and/or  spectrometer  resolution  just  as  In  two  color  wavemixing.  Phase 
matching  Is  easily  achieved  via  a  variety  of  approaches,  e.g.  colllnear, 
closely-packed  three  beam,  combined  planar  BOXCARS  or  combined  planar-folded 
BOXCARS  schemes. 

Due  to  the  spectral  energy  partitioning  Inherent  in  the  use  of  two 
broadband  sources,  the  dual  broadband  signals  are  somewhat  weaker  than  would 
normally  pertain  In  two-color  wave  mixing.  For  dynamic  measurements,  dual 
broadband  CARS  might  be  restricted  to  elevated  pressure  combustion  applica¬ 
tions,  e.g.  gas  turbines,  diesels,  propellant  burning.  At  elevated  pressure, 
the  signal  decrease  accompanying  dual  broadband  CARS  would  be  compensated  by 
the  nominal  quadratic  scaling  of  the  CARS  signal  with  density.  This  new 
technique  Is  currently  being  Implemented  In  the  UTRC  High  Pressure  CARS 
Facility  and  will  allow  simultaneous  measurements  of  all  major  flame  products, 
CO  and  CO^. 
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Dual  broadband  CARS  for  simultaneous,  multiple  species 
measurements 


Alan  C.  Eckbreth  and  Torger  J.  Anderson 


An  approach  to  CARS  is  described  and  demonstrated  which  permits  CARS  to  be  generated  from  a  multiplic¬ 
ity  of  species  simultaneously.  The  technique  employs  two  independent  broadband  Stokes  sources  in  combi¬ 
nation  with  a  pump  laser.  In  addition  to  the  two  separate  two-color  wave-mixing  processes  between  the 
pump  and  Stokes  lasers,  spectrally  resolved  CARS  is  produced  in  a  three-color  process  from  species  whose 
Raman  resonances  correspond  to  the  frequency  differences  of  the  two  broadband  sources.  CARS  is  thus 
derived  from  a  large  number  of  species  simultaneously  removing  the  nominal  limitation  of  CARS  to  interro¬ 
gate  only  one  constituent  at  a  time. 


I.  Introduction 

CARS  (coherent  anti-Stokes  Raman  spectroscopy) 
is  by  now  a  well-established  technique  for  the  remote, 
spatially  and  temporally  resolved  probing  of  instru- 
mentally  hostile  environments  typical  of  practical 
combustion  and  plasma  devices.'"*  A  major  limitation 
of  CARS,  in  addition  to  increased  complexity,  when 
compared  to  spontaneous  Raman  scattering,  is  the  in¬ 
ability  to  interrogate  more  than  one  molecular  species 
at  a  time  as  usually  implemented.  Different  constitu¬ 
ents  are  generally  addressed  sequentially  by  switching 
dye  cells  within  the  Stokes  oscillator  cavity*  or  changing 
the  dye  solution.  In  the  fortuitous  instances  of  closely 
spaced  Raman  resonances,  e.g.,  N2  and  N2O,®  CO2  and 
(>2,  a  single  broadband  dye  laser  allows  more  than  a 
single  constituent  to  be  measured  at  a  time.  Several 
constituents  can  also  be  measured  simultaneously  by 
introducing  a  separate  broadband  Stokes  source  for 
each  species  to  be  measured.  Folded  BOXCARS 
phase-matching  schemes^  ’'  expedite  such  approaches 
and  have  been  employed  with  two  Stokes  sources  to 
monitor  O2  and  N2  simultaneously.®  To  employ  more 
than  two  Stokes  sources  is  clearly  cumbersome  and 
complex.  In  this  paper  we  report  on  a  dual  broadband 
Stokes  approach  to  CARS  which  permits  several,  i.e., 
more  than  just  two,  constituents  to  be  measured  si¬ 
multaneously. 
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In  the  next  section  of  the  paper,  the  basic  dual 
broadband  CARS  concept  is  outlined.  Phase  matching 
all  the  simultaneous  two-  and  three-color  processes  is 
then  described  in  some  detail.  The  paper  concludes 
with  a  discussion  concerning  spectral  positioning  of  the 
Stokes  sources  and  an  experimental  demonstration  of 
simultaneous  CARS  generation  from  postflame  N2, 
CO2,  and  H2O. 

H.  Dual  Broadband  Stokes  Approach 

The  concept,  as  specifically  implemented  here,  is  il¬ 
lustrated  in  Fig.  1.  A  narrowband  pump  laser  (Au]  ;$ 
1  cm"')  wave  mixes  with  two  separate  broadband 
Stokes  sources,  uz  and  (i)'2>  where  u)z  >  ^2.  For  this 
discussion  and  without  loss  of  generality,  the  pump  laser 
is  the  widely  employed,  frequency-doubled  Nd.YAG 
laser  at  532  nm.  Assuming  for  the  moment  that  all 
phase-matching  requirements  are  satisfied,  wi  wave 
mixes  with  uz  and  a>2  each  separately  in  normal  two- 
color  CARS  processes  producing  CARS  at  2«i  -  wz  and 
2wi  —  w'z  for  Raman  resonances  corresponding  to  wi  — 
«2  and  wj  —  <*>2.  respectively.  In  addition,  there  is  the 
three-color  CARS  process,  again  assuming  phase 
matching  is  satisfied,  which  produces  CARS  at  +  uz 
—  w'z  for  Raman  resonances  at  uz  —  u'z-  Depending  on 
the  bandwidths  of  the  broadband  sources,  the  frequency 
differences  uz  ~  (^2  can  span  a  spectral  region  of  many 
hundreds  of  cm"'  permitting  a  multiplicity  of  constit¬ 
uents  to  be  addressed.  Thus  CARS  is  generated  from 
any  constituent  whose  Ramw  resonances  coincide  with 
«i  -  fa)2,  -  u)z,  or  Uz  -  u'z.  Besides  the  three-color 

process  at  wi  +  0)2  -  u'z  for  Raman  resonances  at  uz  - 
u'z,  i.e.,  the  dual  broadband  process,  there  will  also  be 
three-color  wave  mixing  at  the  same  frequency  from 
Raman  transitions  at  wi  -  u'z-  However  for  the  reason 
detailed  below,  the  latter  radiation  will  not  have  high 
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Fig.  1.  Dual  broadband  Stokes  approach  to  simultaneous,  multiple 
species  CARS  measurements. 


resolution  but  will  be  spectrally  smeared.  As  such  it 
could  contribute  to  the  nonresonant  background  and 
complicate  species  concentration  measurements;  thus, 
its  presence  and  importance  need  to  be  scrutinized.  For 
thermometry,  its  effect  could  be  minimized  by 
employing  data  reduction  approaches  which  employ  the 
nonresonant  background  as  a  fitting  parameter.^  In  the 
experiments  to  be  described  later,  its  presence  was  not 
readily  apparent. 

For  the  two,  two-color  cases  and  the  one,  three-color 
process  with  dual  broadband  excitation,  the  instru¬ 
mental  spectral  resolution  is  determined  by  the  pump 
laser  linewidth  and/or  slit  function  even  for  the  three- 
color  mixing.  In  the  latter,  despite  the  employment  of 
the  two  broadband  sources  in  the  three-color  process, 
high  spectral  resolution  results,  because,  as  pointed  out 
by  Yuratich,'®  the  broadband  sources  serve  only  to  ex¬ 
cite  the  vibrational-rotational  Raman  resonances  cor¬ 
responding  to  their  frequency  differences.  Since  the 
Raman  frequencies  are  well  defined  by  the  constituent 
molecules,  spectral  smearing  does  not  occur  even  with 
dual  broadband  excitation.  The  relatively  narrowband 
pump  laser  then  scatters  off  the  resultant  coherent 
Raman  excitation  in  the  medium  resulting  in  CARS  at 
a  spectral  resolution  characteristic  of  the  narrowband 
source  and/or  spectrometer  employed.  On  the  other 
hand,  the  three-color  CARS  process  for  Raman  reso¬ 
nances  at  —  C02  will  result  from  the  broadband  012 
source  scattering  from  the  coherent  excitation  of  those 
modes  and  will  be  spectrally  diffuse.  That  the  three- 
color  CARS  spectrum  with  dual  broadband  excitation 
of  the  Raman  modes  is  indeed  sharp  is  experimentally 
demonstrated  in  Fig.  2  where  the  CARS  spectrum  of 
ambient  N2  is  displayed.  The  experimental  apparatus 
employed  here  is  similar  to  and  is  a  laboratory  version 
of  the  CARS  instrument  described  in  Ref.  4.  The  main 
difference  is  the  employment  of  portions  of  the  532-nm 
beam  to  pump  two  separate  oscillator  cavities  rather 
than  an  oscillator/amplifier  combination  as  commonly 
done.  The  N2  signature  was  generated  in  a  crossed- 
beam  three-color  CARS  process  in  which  two  broad¬ 
band  sources  (150-200-cm“*  FWHH),  centered  at 
~17471  cm~*  (0)2)  and  ~15140  cm“*  («i),  excited  the  N2 
vibrational  fundamental  Q-branch  Raman  resonances 
from  which  the  18,797-cm“’  pump  scattered.  A 
300-mm  focal  length  lens  was  employed  with  the  fol¬ 
lowing  beam  angles  relative  to  the  central  lens  axis;  W], 
1.05®,  0)2, 1.13®  and  0)2, 1.28®.  The  moderate  spectral 
resolution  displayed  (~2.7  cm“*)  is  the  result  of  in¬ 
serting  the  50-M>n  diam  optical  fiber  from  the  CARS 
receiver  directly  into  the  0.75-m  spectrograph  fitted 
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Fig.  2.  Dual  broadband  CARS  signature  of  Na  in  room  air  generated 
from  tbe  three-color  wave  mixing  of  ui],  ua.  and  uia.  Spectral  dis¬ 
persion  is  0.53  cm'Vpixel.  The  spectral  resolution  is  2.7  cm"'.  The 
FWHH  of  the  Na  CARS  spectrum  is  4.5  cm"'. 


with  a  self-scanned  optical  multichannel  detector  (PAR 
model  1420). 

III.  Phase-Matching  Schemes 

To  this  point,  all  the  various  frequency  mixing  com¬ 
binations  have  been  assumed  to  be  phase  matched. 
This  is  readily  achieved  via  a  variety  of  approaches.  In 
gases,  neglecting  dispersion  which  is  generally  valid,  all 
the  various  wave-mixing  combinations  can  be  colli- 
nearly  phase  matched  employing  appropriate  dichroic 
elements  for  beam  juxtaposition  and  signal  separation. 
For  applications  where  high  spatial  resolution  is  re¬ 
quired,  phase  matching  is  possible  with  a  variety  of 
BOXCARS  schemes,®  ’  all  planar,  all  folded,  or  com¬ 
bined  planar-folded.  For  the  initial  experimental 
demonstrations  presented  here,  a  combination  of  planar 
and  folded  approaches  has  been  employed  as  shown  in 
Fig.  3.  Figure  3  depicts  the  laser  beam  arrangement 
incident  on  the  focusing  field  lens  (INPUT)  and  the 
recollimating  field  lens  (OUTPUT)  viewed  against  the 
direction  of  beam  propagation.  U2  is  phase  matched 
with  the  two  wi  beams  via  standard  planar  BOXCARS, 
<1)2  with  «i  by  folded  BOXCARS,  and  the  dual  broad¬ 
band  three-color  process  by  folded  BOXCARS.  In  this 
arrangement,  if  the  individual  two-color  processes  are 
phase  matched,  the  dual  broadband  process  is  phase 
matched  as  well.  This  can  be  shown  by  considering  the 
three  phase-matching  diagrams  shown  in  Fig.  4.  The 
phase-matching  equations  for  the  various  wave-mixing 
combinations  are 


<«'i,a)2 


«»)l,<»2,fa)2 


2ki  cosa  “  *2  cosh  -f  kc  cosc. 

(la) 

*2  8in6  *  hr  sine; 

(lb) 

2^1  cosa  “  *2  oosf  +  kh  cosh, 

(2a) 

k'2  sint  •  kk  sink; 

(2b) 

fe]  COSA  -h  ^2  COsB  “  ^2  COsT  +  kn  cosN, 

(3a) 

k\  sinA  •  k2Binfi, 

(3b) 

h'2  sinT  “  k„  sinN. 

(3c) 
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Fig.  3.  Laser  beam  arrangement  for  the  dual  broadband  CARS  ex¬ 
perimental  demonstrations.  The  beam  geometry  is  shown  just  prior 
to  the  focusing  field  lens  (INPUT)  and  after  the  recollimating  field 
lens  (OUTPl  IT)  viewed  against  the  direction  of  propagation,  wi  and 
0)3  are  mixed  via  planar  BOXCARS,  while  o)i,o)3  and  o)i,o)3.o)2  are  3-D 
phase  matched  (folded  BOXCARS). 


The  photon  energy  conservation  equations  may  be 
written,  assuming  no  dispersion,  as 

2*1  «  *2 -F  Ac.  (4a) 

2ki~kj  +  kh,  (4b) 

*l-F*2**2  +  *n  (4c) 

In  the  small  angle  approximation,  i.e.,  sin9  -►  d,  cos6 
-►  1  -  ff^/2,  Eqs.  (l)-(4)  may  be  recast,  after  some  ma¬ 
nipulation,  to 


6c,  c“  — 6; 


or^  «  th ,  • 


AB~NT,  kiA^h^B,  ^A^~^T^. 

«2 


(Sa.b.c) 


(6a,b,c) 


(7a,b,c) 


2  COLOR  o.’^ .  u)^  2 


3  COLOR,  DUAL  BROADBAND  CARS 

o)^ ,  0)2 


The  small  angle  approximations  are  valid  for  angles  at 


Fig.  4.  Phase-matching  diagrams  for  the  individual  wave  mixing 
combinations  comprising  dual  broadband  CARS. 


least  up  to  10®;  BOXCARS  experiments  very  seldom 
exceed  these  values. 

If  the  CARS  and  Stokes  wave  vectors  are  expressed 
in  terms  of  the  pump  wave  vector  and  Raman  shift  A, 
namely, 

*c  “  *2  “  ~ 

*/■  “  ^1  +  *2  “  ^1  ~  18) 

*,  «  *1  +  A„  »  *1  +  A*  -  Ac 
it  is  easily  shown  for  that 

*f  -  *2*c,  (9a) 

*S-*i*s,  (9b) 

kihi-k'jkn.  (9c) 

To  prove  that  the  three-color  dual  broadband  process 
is  phase  matched  when  the  individual  two-color  pro¬ 
cesses  are  each  separately  phase  matched,  one  proceeds 
as  follows.  Beginning  with  A  ^  a,  one  can  show  that  B 
*  6  by  squaring  Eq.  (7b),  replacing  k\ by  Eq.  (9a),  and 
comparing  the  resultant  expression  with  Eq.  (5b). 
Similarly  T  »  t  can  be  shown  by  substituting  sequen- 
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Fig.  5.  Laser  beam  arrangement  for  dual  broadband  CARS  in  which 
all  wave-mixing  combinations  are  3-0  phase  matched  and  each  a>i 
component  participates  in  the  three-color  process. 


tially  Eqs.  (9c)  and  (9b)  in  Eq.  (7c)  and  comparing  with 
Eq.  (6b).  Thus,  the  three-color  dual  broadband  CARS 
phase-matching  angles  equal  the  two-color  phase¬ 
matching  angles  for  this  combination  of  planar  and 
folded  BOXCARS.  For  very  broadband  sources,  there 
will  be  some  phase  mismatch  for  frequencies  far  from 
line  center.  This  could  result  in  a  signal  decrease  due 
to  lack  of  perfect  phase  matching,  but,  for  small  angle 
beam  crossings,  probably  does  not  result  in  any  spectral 
distortion  within  a  specific  Raman  resonance  band  re¬ 
gion. 

There  are  two  drawbacks  to  the  phase-matching 
scheme  shown  in  Fig.  3.  First,  only  one  of  the  W]  com¬ 
ponents  participates  in  the  three-color  dual  broadband 
process.  Second,  for  high  spatial  resolution  applica¬ 
tions,  the  collinear  CARS  component  generated  in  the 
planar  BOXCARS  process  (wi,a)2)  needs  to  be  elimi¬ 
nated  by  spatial  filtering.  In  Fig.  5,  what  may  well  be 
the  optimal  phase-matching  scheme  for  dual  broadband 
CARS  is  suggested.  In  this  arrangement,  all  the  major 
wave-mixing  processes  are  phase  matched  by  folded 
BOXCARS.  Furthermore,  each  wi  component  par- 


Fig.  6.  Phase-matching  geometry  for  the  laser  beam  arrangement 
of  Fig.  5.  Two-color  processes  are  phase  matched  about  the  central 
lens  axis.  The  three-color  processes  are  phase  matched  about  the 
dashed  axis  residing  in  a  plane  inclined  to  the  uii  beam  crossing 
plane. 

ticipates  separately  in  the  dual  broadband  CARS  pro¬ 
cess,  i.e.,  each  wi  component  scatters  off  the  coherent 
Raman  excitation  established  by  the  nearly  collinear 
W2,(^2  components.  There  wiU  also  be  a  weak,  spectrally 
smeared  collinear  CARS  signal  at  2a)2  -  W2  residing 
along  the  102,0)2  beams.  This  is  easily  eliminated  by- 
trapping  and  is  located  spectrally  at  a  smaller  anti- 
Stokes  shift  from  o)i  than  the  CARS  signal  from  the 
o)i,o)2  interaction. 

It  is  also  possible  to  prove  here,  as  before,  that  the 
three-color  dual  broadband  CARS  process  is  phase 
matched  when  the  individual  two-color  CARS  processes 
are  phase  matched.  The  proof  in  this  instance  is  more 
complicated  since  the  three-color  CARS  process  is  no 
longer  phase  matched  about  the  central  lens  axis  as  in 
the  planar-folded  orientation  of  Fig.  3.  The  phase¬ 
matching  axes  reside  in  planes  tilted  relative  to  the  lens 
axis  as  shown  in  Fig.  6.  For  this  reason,  the  three-color 
processes  will  have  a  slightly  different  interaction 
length,  and  thus  spatial  resolution,  than  the  individual 
two-color  combinations.  In  Fig.  6,  only  the  three-color 
process  for  one  of  the  wi  pump  components  is  shown; 
the  other  process  possesses  an  identical  geometry.  The 
location  of  the  dotted  phase-matching  axis  is  selected 
so  that  angles  A  and  fi  satisfy  the  three-color  phase¬ 
matching  relation  in  Eq.  (7b).  The  proof  proceeds  by 
demonstrating  that  the  angle  T,  solved  for  geometri¬ 
cally,  satisfies  phase-matching  closure  indicated  by  Eq. 
(7c).  We  assume  a  field  lens  of  focal  length  f  much 
larger  than  the  displacements  of  the  laser  beams  from 
the  central  lens  axis  about  which  the  two-color  wave¬ 
mixing  processes  are  phase  matched.  The  displace¬ 
ments  of  the  beams  from  the  central  axis  are  denoted 
by  r,  and  for  r,  «  f,  generally  the  case,  one  may  express 
the  beam  displacements  in  terms  of  the  two-color 
phase-matching  angles  thusly, 

ri  “  erf,  rj  *  b/,  «  tf,  (10a, b,c) 

where  the  subscripts  denote  the  laser  frequencies.  The 
displacements  of  the  beams  from  the  three-color 
phase-matching  axis  are  denoted  by  A,  and  are  related 
to  the  three-color  phase-matching  angles  by 
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Ri-A/.  R2“B/.  R2=Tf.  (11a, b.c) 


One  should  note  that  the  three-color  process  is  not 
planar  but  slightly  folded.  The  three-color  phase- 
matching  axis  is  selected  to  satisfy  the  Eq.  (7b)  phase- 
matcbing  condition.  Thus  the  R2,Ri  partitioning  is 
given  by 


where 


A  ^2 
/?2  B  ki 


(12) 


Ri^R2  =  y/r{Vrl  (13) 

Applying  the  law  of  cosines  to  the  slender,  obtuse 
triangle  formed  by  R2,  R2  and  (rj  —  r2)  one  has 


Ri  =  fi j  -(•  (rj  —  r2)^  -  2/?2('’2  ~  ''2)  cos(180  —  P),  (14) 

where 


cosP  = 


—Ll— 

\'  r'i  ■+  r'i 


(15) 


Dividing  Eq.  (14)  by  P  and  substituting  from  Eqs.  (10) 
and  (11)  transform  Eq.  (14)  to 


T^  =  B^+(t  -  6)2  + 


2B(t  -  6)6 
A  +  B 


(16) 


Dividing  Eq.  (13)  by  /  leads  to 


(A  +  B)2  =  a2  +  62, 


(17) 


or  to 


a 


2  ss 


{k\  +  62)2 


(18) 


after  eliminating  B  and  b  using  Eqs.  (5b)  and  (7b).  In 
Eq.  (16)  B  may  be  expressed  in  terms  of  A,  t  and  b  in 
terms  of  a,  and,  thus,  of  A  leading  to 


species  of  measurement  interest.  For  the  purposes  of 
the  following  discussion,  note  that  most  2  X  Nd;YAG- 
pumped,  broadband  dye  lasers  possess  spectral  profiles 
with  FWHH  of  the  order  of  150-200  cm“'  and  base- 
widths  approximately  twice  this  value.  Binary  dye 
combinations,  which  we  have  not  yet  investigated  ex¬ 
tensively,  can  be  found  which  yield  even  larger  lasing 
bandwidths.  For  combustion  applications,  some  a;2,a;2 
combinations  come  readily  to  mind.  For  0)2,  one  would 
certainly  like  to  exploit  the  broad  frequency  coverage 
provided  by  the  dye  DCM.”  Depending  on  solvent  and 
concentration,  this  dye,  when  2  X  Nd:YAG  numped, 
iases  over  many  hundreds  of  cm“’  and  can  be  used  to 
generate  CARS  from  H2O  and  the  C-H  stretch  modes 
of  a  variety  of  hydrocarbons.  With  DCM  dissolved  in 
DMSO,  we  have  experimentally  seen  Stokes  FWHH  of 
the  order  of  350  cm”'.  Assume  a  base  spectral  coverage 
over  the  3000-3700-cm”’  range.  One  choice  for  W2  then 
is  the  1200-1600-cm~’  region  permitting  CARS  directly 
from  CO2  and  O2.  Dual  broadband  CARS  is  then  ob¬ 
tained  from  species  with  resonances  covered  by  0)2  —  a>2 
or  the  1400-2500-cm~'  range,  e.g.,  N2,  NO,  CO.  N2O,  C2. 
etc.  In  the  flame  work  to  be  presented  subsequently, 
we  centered  0^2  at  H2O  (3657  cm”')  and  0)2  at  1326  cm”' 
so  that  maximum  dual  broadband  CARS  would  be 
generated  from  N2.  Another  choice  for  0)2  would  be  the 
2000-2400-cm”'  region  to  generate  CARS  directly  from 
N2,  CO,  etc.,  with  —  wi  covering  resonances  between 
600  and  17()0  cm”',  e.g.,  CO2,  O2.  If  ^>2  is  set  to  lase 
between  3600  and  4200  cm”'  corresponding  to  H2O  and 
H2  and  a>2  to  13(X)-16(X)  cm”'  (02,C02),  dud  broadband 
CARS  would  be  generated  from  resonances  in  the 
2000-3000-cm”'  range. 

Due  to  the  spectral  energy  partitioning  inherent  in 
the  use  of  broadband  dye  sources  in  wave-mixing  ex- 


fj_  _  +  *1  (1  _ 

162  262  kj 


{ki  +  kt)A^. 


(19) 


Numerical  calculations  indicate,  fortunately,  that  the 
second  term  on  the  right-hand  side  of  Eq.  (19)  is  quite 
small  and  can  be  neglected;  this  is  equivalent  to  drop¬ 
ping  the  second  term  on  the  right-hand  side  of  Eq.  (14). 
After  considerable  manipulation,  using  Eqs.  (9),  Eq. 
(19)  simplifies  to 

62621’^ +  -42.  (20) 

For  typical  values  of  kc  (20, 123  cm”')  and  k„  (21,128 
cm”'),  the  term  in  brackets  is,  from  Eq.  (8),  equal  to 
to  within  0.5%.  Thus  Eq.  (20)  is  equivalent  to 

Mil'S -fe,fe„A2,  (21) 

demonstrating  that  T  does  satisfy  the  phase-matching 
condition  of  Eq.  (7c).  Exact  numerical  calculations 
show  that  the  three-color  process  is  indeed  phase 
matched  indicative  of  the  validity  of  the  approximations 
employed. 

IV.  Dual  Stokes  Spectral  Selection 

The  spectral  positioning  of  the  two  broadband  dye 
laser  sources  clearly  depends  on  the  application  and 


periments,  dual  broadband  CARS  may  be  restricted, 
particularly  for  dynamic  measurements,  to  elevated 
pressure  combustion  applications,  e.g.,  propellant 
burning,  diesel  combustion.  At  elevated  pressures,  the 
signal  decrease  accompanying  the  use  of  dual  broad¬ 
band  sources  is  compensated  by  the  nominal  quadratic 
scaling  of  |  1  ^  with  density.'^  The  exact  scaling  de¬ 

pends  on  the  (detailed  spectroscopy  of  the  particular 
species  in  mind,  e.g.,  transition  spacing,  Raman  line- 
width,  collisional  narrowing  effects,^  etc.  The  signal 
for  any  particular  species  will  also  depend  on  the  (02  ~ 
u>2  overlap  integral  with  Uie  Raman  resonances  excited. 
For  example,  if  o>2  and  oiz  have  the  same  bandwidth  and 
the  center  frequencies  U2(0)  and  w^iO)  coincide  with  a 
particular  Raman  resonance,  all  the  energy  in  each  laser 
is  available  to  excite  those  particular  Raman  resonances. 
This  is  unlike  normal  broadband  CARS  with  a  nar¬ 
rowband  pump  source  in  which  only  a  fraction  of  the 
bro;idband  dye  laser  energy  participates  in  the  resonant 
wave  mixing.  Thus,  the  signal  decrease  with  dual 
broadband  CARS  relative  to  normal  two-color  CARS 
is  probably  not  as  large  as  one  might  originally  antici- 
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Fig.  7.  Simultaneously  generated  dual  broadband  CARS  signatures 
from  Ns.  COj,  and  HsO  in  the  pdetflame  zone  of  an  "IVOO  K  premixed 
CH4-air  flame.  The  dispersions  are  ~0.53  cm'Vpixel  for  Ns,  0.59 
cm'Vpixel  for  HsO,  and  0.47  cm'Vpixel  for  CO3.  In  CO2  both  the 
vi  and  2vs  bands  are  displayed. 

pate.  In  our  initial  experiments,  for  example,  with 
W2(0)  -  wsiO)  =  An2,  dual  broadband  CARS  from  Na  in 
the  postflame  zone  of  an  atmospheric  pressure,  pre¬ 
mixed  CH4-air  flame,  was  readily  recorded  with  time 
averaging  over  600  laser  pulses  or  30  sec  and  a  spatial 
resolution  of  several  millimeters  as  shown  in  Fig.  7. 
Also  shown  there  are  the  postflame  CARS  spectra  from 
CO2  and  H2O  simultaneously  generated  with  the  N2 
CARS  signal  and  sequentially  recorded  by  merely 
tuning  the  spectrograph.  These  of  course  could  be  si¬ 
multaneously  recorded  employing  a  spectrograph  with 
less  dispersion  than  used  here  or,  preferably,  placing  the 
spectral  images  at  high  dispersion  on  different  portions 
of  either  a  linear  or  2-D  array  detector  using  various 
alignment  schemes.  In  practical  measurement  appli¬ 
cations,  many  of  the  simultaneously  generated  CARS 
signals  would  be  recorded  as  spectrally  integrated 
photomultiplier  outputs. 

When  the  beam  crossing  angles  are  small,  all  the 
CARS  beams  generated  can  be  conveniently  focused 
into  the  spectrograph  or  receiver  optical  fiber^  with  a 
single  lens.  However,  as  the  beam  crossing  angles  in¬ 
crease,  the  spatial  separation  of  the  CARS  beams  may 
lead  to  incidence  angles  in  excess  of  that  acceptable 
based  on  the  //No.  of  the  spectrograph  or  numerical 
aperture  of  the  fiber.  In  those  instances,  separate 
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capture  of  the  individual  CARS  beams  can  be  employed 
with  subsequent  optical  repositioning  to  meet  the  //No. 
constraints  of  the  optical  detection  system. 

It  is  also  interesting  to  speculate  on  the  single  pulse 
quality  of  dual  broadband  CARS  spectra  when  detector 
shot  noise  is  not  problematical.  Normal  two-color 
single  pulse  CARS  spectra  are  affected  by  the  amplitude 
ripple  on  the  broadband  Stokes  dye  laser.^  This  leads 
to  distorted  single  pulse  CARS  spectra  and  a  decrease 
in  the  accuracy  of  single  shot  measurements.  For  two 
statistically  independent  dye  sources,  where  many 
frequency  combinations  drive  each  individual  Raman 
resonance,  one  might  expect  dual  broadband  CARS 
single  pulse  quality  to  be  considerably  better.  This 
would  arise  from  an  averaging  in  effect  over  the  random 
amplitude  profiles  of  the  two  broadband  sources  which 
should  lead  to  smoother  single  pulse  spectra.  This  as¬ 
pect,  as  well  as  others,  will  be  reported  as  investigations 
of  dual  broadband  CARS  progress. 

V.  Summary 

We  have  described  a  readily  implemented,  dual 
broadband  approach  to  CARS  which  permits  many 
molecular  species  to  be  monitored  simultaneously  thus 
eliminating  one  of  the  major  drawbacks  to  CARS  di¬ 
agnostics.  The  precise  range  of  sensitivity  and  appli¬ 
cability  of  the  technique  will  become  clearer  with  fur¬ 
ther  investigations  in  this  regard  and  engineering  de¬ 
velopments. 
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